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Q _ Abstract 

We construct spectral triples for path spaces of stationary Bratteli diagrams and 
£ — ■ study their associated mathematical objects, in particular their zeta function, their 

CN1 . heat kernel expansion and their Dirichlet forms. One of the main difficulties to 

properly define a Dirichlet form concerns its domain. We address this question in 
particular in the context of Pisot substitution tiling spaces for which we find two 
types of Dirichlet forms: one of transversal type, and one of longitudinal type. Here 
the eigenfunctions under the translation action can serve as a good core for a non- 
' trivial Dirichlet form. We find that the infinitesimal generators can be interpreted 

as elliptic differential operators on the maximal equicontinuous factor of the tiling 
dynamical system. 
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1 Introduction 

Even though noncommutative geometry [9| was invented to describe (virtual) noncom- 
mutative spaces it turned out also to provide new perspectives on (classical) commutative 
spaces. In particular Connes' idea of spectral triples aiming at a spectral description of 
geometry has generated new concepts, or shed new light on existing ones, for topological 
spaces: dimension spectrum, Seeley type coefficients, spectral state, or Dirichlet forms 
are notions which are derived from the spectral triple and we will talk about here. In- 
deed, we study in this paper certain spectral triples for commutative algebras which are 
associated with stationary Bratteli diagrams, that is, with the space of infinite paths on 
a finite oriented graph. Such Bratteli diagrams occur in systems with self-similarity such 
as the tiling systems defined by substitutions. 

Our work can be seen as part of a largeer scheme to study self-similar structures with 
the help of spectral triples. This has mainly been done for fractals |12l [T3"l [TJ |8] and 
recently also for tilings |20[ 116] . There has also been work on spectral triples for AF- 
algebras [1], which are (mostly non-commutative) algebras that may be very effectively 
described by Bratteli diagrams. Although the algebras we look at here are commutative 
AF-algebras our approach to constructing a spectral triple is quite different from the one 
in [3] and does not easily generalise to non-commutative AF-algebras. 

Our construction derives from earlier onces for metric spaces which go under the 
name "direct sum of point pairs" [5j or "approximating graph" |19] . suitably adapted 
to incorporate the self-similar symmetry. The construction is therefore more rigid. The 
so-called Dirac operator D of the spectral triple will depend on a parameter p which 
is related to the dilation of the self-similarity. We observe a new feature which, we 
believe, ought to be interpreted as a sign of self-similarity: The zeta function is periodic 
with purely imaginary period Correspondingly, what corresponds to the Seeley 

coefficients (in the case of manifolds) in the expansion of the trace of the heat-kernel 
e -tD j g hgj-g given by functions of logt which are j^^-periodic. This has consequences 
for the usual formulae for tensor products of spectral tiples. If we take the tensor product 
of two such triples and compare the spectral states 71,2 for the individual factors with the 
spectral state T of the tensor product, then a formula like T{A\ <g> A2) = T\ (^1)72(^2) 
will not always hold due to resonance phenomena of the involved periodicities. 
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The heat kernel we were referring to above is the kernel of the semigroup generated 
by D 2 and hence does not involve the algebra. But the spectral triple gives in principle 
rise to other semigroups whose generators may be defined by Dirichlet forms of the type 
Q{fi9) = T([D, f]*[D,g]). Here /, g are represented elements of the algebra and T a 
state. We will take for T the spectral state, a common choice, but not the only possible 
one. Pearson-Bellissard |20] . for instance, choose the standard operator trace. There is 
however a difficulty, namely it is a priori not clear what is the right domain of definition 
of Q. As we will conclude from this work, the choice of domain is crucial and needs 
additional ingredients. This is why we can only discuss rigorously Dirichlet forms and 
Laplacians (their infinitesimal generators) in the second part of the paper, when we 
consider our applications. 

Our main application will be to the tiling space of a substitution tiling. In this case 
the finite oriented graph defining the spectral triple is the substitution graph. Moreover, 
the spectral triple is essentially described by the tensor product of two spectral triples 
of the above type, one for the transversal and one for the longitudinal direction. There 
will be thus two parameters, pt r and p\ g . The additional ingredients which will allow 
us to define a domain for the Dirichlet form Q are the dynamical eigenfunctions of the 
translation action on the tiling space. We have to suppose that these span the Hilbert 
space of L 2 -functions on the tiling space and we are thus lead, by Solomyak's theorem 
|21| . to consider Pisot substitutions. This means that the Perron- Frobenius eigenvalue of 
the substitution matrix is 9 d , where 9 is a Pisot number and d the dimension of the tiling. 
Our main result about the Dirichlet form can then be qualitatively explained as follows: 
In order to have non-trivial Dirichlet forms the parameters pt r and pi g have to be fixed 
such that pig = 9~ 1 and pt r = \9'\ where 9' is an algebraic conjugate to 6, distinct from 9, 
with maximal modulus. The modulus \9'\ is strictly smaller than 1 by the Pisot-property 
and larger than 9~ l . It is equal precisely for quadratic unimodular Pisot numbers. It 
then follows that the Laplacian defined by the Dirichlet form can be interpreted as an 
elliptic operator on the maximal equicontinuous factor of the translation action on the 
tiling space. 

Summary of results After a quick introduction to spectral triples we are first con- 
cerned with the properies of their zeta functions in the case that the expansion of the 
trace of the heat kernel e~ tD is not simply an expansion into powers of t but of the type 

Tr(e-' z?2 )^°/(-logt)t a , (1) 

with K(a) < 0, and / : M + — > M a bounded locally integrable function such that 
lim s _>.o+ s £[fa](s) exists and is non zero, where £ is the Laplace transform. A non- 
constant / in that expansion has consequences which we did not expect at first. We are 
lead in Section[2l!]to study classes of operators on B(Sj) which have a compatible behav- 
ior. An operator A € B(F>) is weakly regular if there exists a bounded locally integrable 
function Ja '■ — > M, for which lim s _;.Q+ s£[fA](s) exists and is non zero, such that 

Tr(e- iD2 ^)^°/A(-logt)^, (2) 
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where a is the same as in equation ([TJ. For such operators, the spectral state does not 
depend on a choice of a Dixmier trace and is given by: 

7*4) - Hm^M, 
s~^o+ C[f\{s) 

where / is the same as in equation ([1]) , see Lemma 12.31 We also define strongly regular 
operators, for which one has in particular Ja = T(A)f in equation ([2]) (see Lemma l2,4p . 
Regular operators have an interesting behavior under tensor product, which we will use 
in the applications to tilings. If the spectral triple is a tensor product: (A,S),D) = 
(A\ ®A 2 ,S)i ®Sj 2 ,Di <g) 1 + D2), where x is a grading on (Ai,Sji,Di), then one has: 

T(A 1 ® A 2 ) = hm , 
s ^ + £[/i/ 2 J(s) 

where /, is as in JI]) for Di , and as in ([2]) for , for each factor i = 1, 2 of the tensor 
product, see Lemma 12.71 In general, only if both A\ and j4 2 are strongly regular for 
the individual spectral triples the state will factorize as T(A\ ® A 2 ) = 7% (A±) 72^4.2 )■ 
Here % denotes the spectral state of (Ai,$)i, Di), i = 1,2 (see Corollary 12. 8|) . It is easy 
to build examples for which this equality fails for more general operators: for example 
Ti(Ai) = T 2 (A 2 ) = and T{A 1 ® A 2 ) ^ 0. 

In Section [3] we study spectral triples associated with a stationary Bratteli diagram, 
that is for the C* -algebra of continuous functions on the Cantor set of (half-) infinite paths 
on a finite oriented graph. These depend on the matrix A encoding the edges between 
two levels in the diagram (called here a graph matrix and assumed to be primitive), a 
parameter p E (0, 1) to account for self-similar scaling, and a horizontal structure % 
(a set of edges linking the edges of the Bratteli diagram). We determine the spectral 
information of such spectral triples. In Theorem 13.31 we derive the Connes-distance, 
and show under which conditions it yields the Cantor topology on the path space. We 
compute the zeta-function £(z) = Tr (\D\~ Z ) and the expansion of the heat-kernel. 
Theorem (Theorems 13.41 and !3.6l in the main text.) Consider a spectral triple associated 
with a stationary Bratteli diagram with graph matrix A and parameter < p < 1. Assume 
that A is diagonalizable with eigenvalues Ai, ■ • • , A p . 

• The zeta-function £ extends to a meromorphic function on C which is invari- 
ant under translation z >— >■ z + j^-- It has only simple poles and these are at 

l0g - log 2 p" fc » k E Z,j = 1, . . .p. In particular, the spectral dimension (abscissa of 
convergence of X) is equal to sq = '"^g p F , where A PF is the Perron- Frobenius eigen- 
value of A. The residue at the pole log A 7 + 27Ttk ^ q{ ven b v c n x i 

J f — log p z a _ | g p 

• The Seeley expansion of the heat-kernel is given by 

Tr(e"^ 2 )= £ ^p_ 21ogp , log A,(-logt) tS* +C%-=££- + h(t) 



where h is entire, p r ^ a is an r -periodic smooth function, and jo is such that A 
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The constants are given in ([10p . they depend on the choice of horizontal edges H. 
The function p r ^ a is explicitly given in equations (|15p and (|16p . and its average over a 
period is p rj0 = -T(-). If A is not diagonalizable then £ has poles of higher order and 
the heat-kernel expansion is more involved (with powers of log(t) depending on the order 
of the poles), see Remarks 13.51 and 13.71 

In Section U] we apply our findings to substitution tiling spaces £l a . We consider 
geometric substitutions of the simplest form, as in |14| . which are defined by a decom- 
position rule followed by a rescaling, i.e. each prototile is decomposed into smaller tiles, 
which, when stretched by a common factor 9 > 1 (the dilation factor) are congruent to 
some original tile. The result of the substitution on a tile is called a supertile (and by 
iteration then an n-th order supertile). If one applies only the decomposition rule one 
obtains smaller tiles. These we call microtiles. 

The substitution induces a hyperbolic action on Our approximating graph will 
be invariant under this action. But £l a carries a second action, that by translation of the 
tilings. Although the translation action will not play a direct role in the construction of 
the spectral triple, it will be crucial in Section 14.71 to define a domain for the Dirichlet 
forms. 

The approximating graph for Q, a is constructed with the help of doubly infinite paths 
over the substitution graph. Half infinite paths describe its canonical transversal. We 
use this structure to construct a spectral triple for C(£l a ) essentially as a tensor product 
of two spectral triples, one obtained from the substitution graph and the other from the 
reversed substitution graph. Indeed, the first of the two spectral triples describes the 
transversal and the second the longitudinal part of Q a . Since the graph matrix of the 
reversed graph is the transpose of the original graph matrix we will have to deal with 
only one set of eigenvalues Ai, • • • , X p . It turns out wise, however, to keep two dilation 
parameters pt r and pi g as independent parameters, although they will later be related to 
the dilation factor 9 of the substitution. We obtain: 

Theorem (Theorem 14.71 in the main text.) The spectral triple for C{£l a ) is finitely 
summable with spectral dimension 

d log 9 d log 9 
so = — ; 1 ; , 

- log ptr ~ log pig 

which is the sum of the spectral dimensions of the triples associated with the transversal 
and to the longitudinal part. 

The zeta function £(z) has a simple pole at sq with positive residue. 

The spectral measure is equal to the unique invariant probability measure on Q a . 

We discuss in Section 14.61 the particularities of Pisot substitutions. These are sub- 
stitutions for which the dilation factor 9 is a Pisot number: an algebraic integer greater 
than 1 all of whose Galois conjugates have modulus less than 1. Their dynamical system 
(J7 CT ,R rf ) factors onto an inverse limit of dJ-tori, its maximal equicontinuous factor E, 
where J is the algebraic degree of 9. The substitution induces a hyperbolic homeomor- 
phism on that inverse limit which allows us to split the tangent space at each point into 
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a stable and an unstable subspace, S and U. The latter is (i-dimensional and can be 
identified with the space in which the tiling lives. S can be split further into eigenspaces 
of the hyperbolic map, namely S = Si + S2 where S2 is the direct sum of eigenspaces 
to the Galois conjugates of 9 which are next to leading in modulus, that is have maxi- 
mal modulus among the Galois conjugates which are distinct from 0. This prepares the 
ground for the study of Dirichlet forms and Laplacians on Pisot substitution tiling spaces 
in Section T4.71 The main issue is to find a domain for the bilinear form on C(£l a ) defined 



Q decomposes into two forms Qt r and Qi g , a transversal and a longitudinal one, which 
turn out to be Dirichlet forms on a suitable core, once the parameters have been fixed to 
Ptr = |^2| an d pig = where 62 is a next to leading Galois conjugate of 0. Our main 
theorem is the following: 

Theorem (Theorem 14.121 in the main text.) Consider a Pisot substitution tiling of M. d 
with Pisot number 8 of degree J. Suppose that the tiling dynamical system has purely 
discrete dynamical spectrum. Let 9%,---6l be the subleading conjugates of so that in 
particular 9j = |#2|e* aj for 2 < j < L < J. Assume that for all j 7^ f one has 
ctj — ctji + 2nk + 27r 1 °p|^ k' 7^ , \/k, k' £ Z. Then the space of finite linear combinations 
of dynamical eigenf unctions is a core for Q on which it is closable. Furthermore, Q = 

Qtr + Qig, and Qt r jig has generator ^-tr/lg = J2heH tr / l ° ^tr/lg 9^ ven 071 an eigenfunction 
fg to eigenvalue (3 by 



We explain the notation as far as possible without going into details. The horizontal edges 
encode vectors of translation between two tiles in some supertile. Whereas a transversal 
horizontal edge h € r H} r encodes a vector which should be thought of as a return 
vector between supertiles of a given type sitting in an even larger supertile-and hence as 
large vectors-, a longitudinal horizontal edge h G % l9 stands for a vector of translation 
ah in a microtile between even smaller microtiles-and so are small. t s 2nA is the tile (up to 
similarity) from which the translation encoded by h starts and freq(t s 2^)) its frequency 
in the tiling. By definition an eigenvalue is an element (3 G R rf * for which exists a function 
f/3 '■ tie C satisfying fp(ui + a) = e 2m ^^ fp(ui) (we write the translation action simply 
by (a,u;) 1— > OJ + a) but the geometric construction discussed in Section [4.61 allows us to 
view j3 also as an element of S © U so that /3(a) = (a,/3) ((•, •) is a scalar product on 
S © U). Here we wrote a for the vector in U corresponding to a via the identification of 
U with the space in which the tiling lives. Finally * : U — > S2 C S is the reduced star 
map. This is Moody's star map followed by a projection onto S% along S\. The values 
for the constants Q r and c\ g are given in (|4"T|) and (I38p . 

The transverse Laplacian can therefore be seen as a Laplacian on the maximal 
equicontinuous factor E. The longitudinal Laplacian can be written explicitly on Q a , 
and turns out to be a Laplacian on the leaves, namely it reads Ai g = q g V| JCV i g , where 
V/ g is the longitudinal gradient, and K, a tensor, see equation (|39p . 



by 



Q(f,g) = T{[D,7T(f)}*[D,7T 



(9)}) ■ 




cig{2Tr) 2 {req(t s 2 {h) )f3(a h ) 2 f 8 , 
Ctr{2Tr) 2 freq(t s 2 ih) )(f h *,^) 2 ld jfg. 
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2 Preliminaries for spectral triples 



A spectral triple (A, 5), D) for a unital C*-algebra A is given by a Hilbert space fj carrying 
a faithfull representation tt of A by bounded operators, and an unbounded self-adjoint 
operator D on 9) with compact resolvent, such that the set of a G A for which the 
commutator [D, ir(a)] extends to a bounded operator on S) forms a dense subalgebra 
.4,0 C -4. The operator L> is referred to as the Dirac operator. The spectral triple 
(A,$),D) is termed even if there exists a Z/2-grading operator \ on ■£> which commutes 
with vr(a), o£i, and anticommutes with D. 

We will here consider the case of commutative C*-algebras, A = (C(X), || • |[oo), of 
continuous functions over a compact Hausdorff space X with the sup-norm, so we may 
speak about a spectral triple for the space X. The spaces we consider will be far from 
being manifolds and our main interest lies in the differential structure defined by the 
spectral triple. More specifically we restrict our attention to the Laplace operator(s) 
defined by it. The approach to defining Laplace operators via spectral triples has been 
considered earlier, for fractals by |12| \13[ ?] and for ultrametric Cantor sets and tiling 
spaces in |20l US! US] . 

2.1 Zeta function and heat kernel 

Since the resolvent of D is supposed compact Tr can be expressed as a Dirichlet 

series in terms of the eigenvalues of \D\^\ The spectral triple is called finitely summable 
if the Dirichlet series is summable for some s G R and hence defines a function 

C(z) = Tr(|Dr), 

on some half plane {z£C: > sq} which is called the zeta-function of the spectral 
triple. The smallest possible value for sq in the above (the abscissa of convergence of the 
Dirichlet series) is called the metric dimension of the spectral triple. We call £ simple 
if lim s _ !>s +(s — sq)C(s) exists. This is for instance the case if £ can be meromorphically 
extended and then has a simple pole at sq- We will refer then to the meromorphic 
extension also simply as the zeta function of the triple. 

Another quantity to look at is the heat kernel e~ tD of the square of the Dirac 
operator. Thanks to the Mellin transform 

T(s)fj,- 2s = / e-^t^dt, 
Jo 

where \i > and T(s) = J + °° e~ t t s ~ 1 dt is the Gamma function, one can relate the zeta 
function to the heat kernel as follows: 

r(s)C(2s)= / Tr(e- w2 )t s ~ l dt. 
Jo 

1 For simplicity we suppose (as will be the case in our applications) that ker(D) is trivial, otherwise we 
would have to work with Tr ker (.D)x (|-D|~ S ) or remove the kernel of D by adding a finite rank perturbation. 
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This of course makes sense only if e~ is trace class for all t > 0, which is anyway a 
necessary condition for finite summability. Notice that the trace class condition implies 
also that s i— > j?° Tr (e~ tD )t s ~ 1 dt is holomorphic for all 5 > 0. 

The above formula is particularity usefull if one knows the asymptotic expansion of 
Tr(e~ tD2 ) at t — >■ 0, or only its leading termH It is well known that the form of the 
asymptotic expansion is related to the singularites of the zeta-function [10\ 115] , For 
instance, an expansion of the form 

Tr(e- tD2 ) = Y J C a t a + h(t), 

a 

with K(a) < 0, c a 6 C, and h a function which is bounded at (in particular without 
logarithmic terms like logi) implies that the zeta-function has a simple pole at —2a 
with residue equal to 2c a /T(— a) and is regular at |10j . We will see however, that the 
situation is quite different in our case where we have to replace c a by functions that are 
periodic in log t. We assume therefore in the sequel that the asymptotic behaviour of the 
trace of the heat-kernel is given by 

Tr(e-* D2 )^°/(-logtr (3) 

where a < and / : IR-° — )■ R is a bounded, locally integrable function for which 
lim s _>o sC[f](s) exists and is different from 0, where C[f] is the Laplace transform of /: 

poo 

£[/](*):= / f(x)e-«°dx. (4) 
J o 

Lemma 2.1. If the trace of the heat-kernel satisfies |3P then £ is simple, has abscissa of 
convergence 

sq = —2a and -T(— ) lim (s — so)C( s ) = hms£[/](s). 

2 2 s ^s+ s^o 

If, moreover, C[f](s) admits a meromorphic extension with simple pole atO and Tr (e~ tD2 ) — 
/(— logt)t a = 0(ir) (with f3 > a) then ((s) has a simple pole at sq = —2a and hence 
|r(f )Res(C,s ) = Res(£[/],0) = lim s ^ sC[f](s). 

Proof. We adapt the arguments of |15| . Let h(t) = Tr(e~ tD2 ) — /(— log t)t a and M = 
sup x Then for all e > exists 6 < 1 such that \h(t)\ < eMt a if t < 5. In particular, 

H 5 (s) := £h(t)t s ~ 1 dt satisfies \H s {s)\ < eMf^, provided a + s > 0. Now, again for 
a + s > 

T(s)C(2s)= TT{e- w2 )t s ~ 1 dt= f(-logt)t a+s - 1 dt + H s (s) + g e ( S ) 
Jo Jo 

2 A function / : R >0 — > C is asymptotically equivalent to g : R >0 — > M at t — > 0, written / g, if 
\f-9\ = o(\g\). f = O(g) means that 3M > 0, 35 > 0,V0 < t < <5 : \f(t)\ < M\g(t)\, and / = o(g) means 
that Ve > 0,35 > 0,V0 < t < 8 : \f(t)\ < e\g(t)\. 
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where g € (s) = J 5 °° Tr (e tD2 )t s 1 dt is holomorphic in s. This shows that C(2s) is finite 
for s > —a. Furthermore 

lim (a + s) [ f(-logt)t a+s ~ 1 dt = lim s [ S f(- T )e TS dT = lim sC[f](s) 

where we have used in the last equation that lim s _ > o+ s Jio g< 5 f(~ T ) eTS( ^ T = 0- Since e > 
is arbitrary we conclude that 

lim (a + s)r(s)C(2s) = lim sC[f]{s). 

s— >— a+ s— >0+ 

Hence so = —2a is the abscissa of convergence. 

Now if h(t) is of order t& we can find M > and 5 > such that |/i(t)t _/3 | < M if 
< i < <5. If /3 > a the function f i— >■ is integrable on (0, 5) as long as s lies in 

a sufficiently small neighbourhood of —a. Since s <— > t s+ ^~ 1 is holomorphic for all t > 
we find that Hg(s) = jQ(h(t)t~^)t s+ ^~ 1 dt is holomorphic near s = —a which shows that 
(a + s)£(2s) is holomorphic there, too. Thus £ has a simple pole at —2a and we have 
the above stated formula for its residue. □ 

Remark 2.2. If f(r) = e taT then C[f](s) = Thus if f is the restriction of a 

periodic function of class C 1 then upon using its representation as a Fourier series we 
see that s£.[f](s) extends to an analytic function around and 



lim sC[f] (s) = f 

s— >0 



the mean of f . 



2.2 Spectral state 

Given a bounded operator A on S) such that |L>| -S °^4 is in the Dixmier ideal we consider 
the expression 

T(A)=Tr u (\D\-' a A)/Tr u (\D\-°) 

which depends a priori on the choice of Dixmier trace Tr w . With a little luck, however, 
lim , + 77-rTr exists and then [9| 

T(A)= lim -Llr(|£>rM). 

s->s+ CO) 

We provide here a criterion for that. Note that the Mellin transform allows us to write 

1 r°° 

Tr(\DrA) = —j o Tr{e- tD2 A)ti~ l dt. 

We call A € B(S)) strongly regular if there exists a number ca such that 
Tr (e~ tD2 A) - c A Tr (e~ w2 ) = o(Tr (e~ w2 )) . 
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If ca 7^ one can thus say that Tr (e~ tD A) ~ c^Tr (e~ tD ). In the context in which 
the heat kernel satisfies Q it is useful to consider the notion of weakly regular operators 
A E B(Sj). These are operators which satisfy 

Trie-^A)^ f A (-logt)t a (5) 

where a is the same as in Q and Ja ■ R-° — > C is a bounded, non-zero, locally integrable 
function for which lim^o sC[/a](s) exists. Clearly, strongly regular operators are weakly 
regular and /a = ca/ in this case, where / is given in equation ([3]) (one actually has 
ca = T(A) , see Corollary [27 



Lemma 2.3. Assume that the trace of the heat-kernel satisfies ([TJj and that A € B(Sj) is 
weakly regular, that is, satisfies Then lim s _ >s + -^^Tr (\D\^ S A) exists and is equal to 

T(A) = lim C[fA]{s) 



>o £[f](s) ' 

Proof. Under the hypothesis, for all e we can find a 8 such that, if s > sq, 



,5 



f (Tr(e- tD2 A)-f A (-logt)t a )r2- 1 dt<e[ \f A (-\ogt)\t a+ ^dt<eM A 
Jo Jo 



where Ma is an upper bound for \f A (— log t) \ . Since Tr (e tD2 )tz 1 dt and hence also 
f 5 °° Tr (e- w2 A)t^ ~ l dt are finite for all 5 > we get (a = -f ) 



lim 



/ Tr (e- w2 A)ti~ l dt - \ f A (- log^t^^dt 



4 r(f )C( S ) 

Notice that £ logt)^ 1 ^ = £[/a](s). 

Since e was arbitrary we conclude that 

S -> S +C(s) u 1 J r(f)c( s ) -o £[/]( s ) 

□ 

Corollary 2.4. If A £ B(Sj) is strongly regular, then Tr (e~ iI)2 ^) *~° T(A)Tr (e'^A). 
In other words, the functions in equations ([3]) and §5§ satisfy f a = T(A)f. 

Proof. If A is strongly regular, then it is also weakly regular with f A = CAf ■ The Laplace 
transform is linear so = c^^C [/](«), and Lemma [2.31 then implies ca = T{A). □ 

Order the eigenvalues of \D\ increasingly (without counting muliplicity) and let F n 
be the n-th eigenspace of \D\. 

Corollary 2.5. Let A G B(Sj) and A n = Tr g^^" ) . // the limit 

A = lim A n 

n— >oo 

exists then A is strongly regular and T(A) = A. 
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Proof. Write c n = e dimi^ where [i n is the n-th eigenvalue of \D\. One has 
Tr(e- w2 A)-ATr(e- w2 ) = £(i n - A)c n 



n>l 



Now fix an e > 0, and choose and integer N £ such that \A n — A\ < e for all n > N e . Then 
the series of the r.h.s. can be bound by (sup n \A n — A\) J2 n <N f c ™ + e Sn>7V E c ™- Using 
^ n> jy c n < Tr (e~ <£>2 ) we find that for all e exists C e such that 



Tr{e- w2 A) - ATt {e" tD ) 



< C e + eTr (e 



-tD 2 



Since Tr (e - *^ ) tends to +00 if t tends to this shows that Tr (e"* D2 A) '~° f A (- log t)t a 
with /a = Af. Applying Lemma [2,31 we see that T(A) = A. □ 

In the commutative case, when A = C(X) for a compact Hausdorff space X, we 
are particularly concerned with operators of the form A = 7r(/), for / £ C(X) or for 
any Borel- measurable function f on X. By Riesz representation theorem the functional 
/ h-» T(7r(/)) gives a measure on X, which we call the spectral measure. 



2.3 Laplacians 

It is tempting to define a quadratic form by 

Q(a,b) = T([D,Tr(a)]*[D,n(b)]) (6) 

on elements on which this expression makes sense. It is however not so easy to determine 
a domain for such a form. Our interest here lies in the commutative case, A = C(X), 
and our question is as follows: Let [i be the spectral measure on X and consider the 
Hilbert space 1?{X) = L 2 (X,fi). Notice that this is usually not the Hilbert space of 
the spectral triple. We embed _4o (continuous functions having bounded commutator 
with D) into L 2 (X), assuming that the measure is faithful. Does there exist a core (a 
dense linear subspace of L 2 (X)) which is contained in .4,0 and on which Q is well-defined, 
yielding a symmetric closable quadratic form which satisfies the Markov property? If T> 
is the domain of the closure of Q then the general theory garanties the existence of a 
positive operator A such that Q(f,g) = (f,Ag) and this operator generates a Markov 
process on L 2 (X). We refer to A as a Laplacian. We emphasize that A will depend on 
the choice of a domain. We won't be able to give general answers here, but we look at 
specific examples related to self-similarity. 



2.4 Direct sums 

The direct sum of two spectral triples («4i,i5i,-Di) and (A2,&2, D2) is the spectral triple 
(A,Sj,D) given by 

A = Ai <S)A 2 , i3=iliei32, D = D 1 (BD 2 . 
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with direct sum representation. The zeta function ( of the direct sum is clearly the sum 
of the zeta functions Q of the summands and thus its abscissa of convergence sq is equal 
to the largest abscissa of the two zeta functions Q. Let T denote the spectral state of 
the direct sum triple and % those of the summands and assume that all zeta-functions 
are simple. Then, for regular operators A\ and A 2 we have 

T(A) = (ciT(Ai) + c 2 T(A 2 )) (7) 

Cl + C 2 

where Cj = lim s _ 5-s +(s — So)d(s). Notice that c\ = if the abscissa of convergence of 
(1 is smaller than that of ( 2 , an d similiarily with 1 and 2 exchanged. Notice also that 
ci + c 2 = lim s ^ s +(s - s )CO)- 

It is pretty clear that the quadratic form Q used to define the Laplacian is the sum 
of the quadratic forms of the individual triples, again leaving questions about its domain 
aside. 



2.5 Tensor products 

The tensor product of an even spectral triple (Ai,S)i, D\) with grading operator x an d 
another spectral triple (A 2 , $) 2 , D 2 ) is the spectral triple (A,f),D) given by 

A = Ai®A 2 , i}=i0i®f3 2 , D = Di <g>l + x®D 2 . 

Notice that D 2 = D\ ® 1 + 1 (g) Djj . It follows that the trace of the heat kernel is multi- 
plicative: Tr(e" tD ) = Tr S j 1 (e- tD i)Tr S j 2 (e- tD 2). This allows one to obtain information 
on the zeta function, the spectral state, and the quadratic form Q of the tensor product. 

Lemma 2.6. Suppose that Tr ^ (e~ tD i ) and Tr^ 2 (e~ tD 2) satisfy ([3jj with f = fi and f = 
f 2 , respectively. Suppose that lim s ^Q sC[f if 2 ](s) exists and is non-zero. Then the metric 
dimension sq of the tensor product spectral triple is the sum of the metric dimensions of 
the factors, and the zeta function £ of the tensor product is simple with 

±rA hm . (s - S o)CW = lim^[/i/ 2 ](a). 

Proof. Due to the multiplicativity of the trace of the heat kernel we have 

Tr (e- w2 ) *2? h{- log t)f 2 {- log t)t a ^ 
and hence the result follows from Lemma [2JJ □ 



Lemma 2.7. Assume the conditions of Lemma \2.(A Let A\ € B{S)\) and A 2 € B(S} 2 ) be 
weakly regular with functions f^ and f^_ 2 . Then 

T(A 1 ®A 2 )= hm C ^ A ^ A ^ 



s^o+ C[hf 2 ](s) 
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Proof. Let e > and choose 5 such that |Tr ^. (e tD i Ai) — f Ai (— log t)t Qi | < et ai provided 
< t < 5. Then 

S 5 

[ Tt ?) {e- tD2 A l ®A 2 )t s - 1 dt = [ Tr^ 1 (e- tD iA 1 )Tr si2 (e- w 2A 2 )t s - 1 dt 
Jo Jo 

= [ (/ Al (-logt)/ A2 (-logt) + 0(e))r i+a2+s - 1 di 
Jo 

= C[f Al f A2 ](a 1 + a 2 + s) + 0(e) 



ati + a 2 + s 

from which the result follows by arguments similar to the above. □ 
Corollary 2.8. Let A\ and A 2 be weakly regular operators. 

(i) If A\ is strongly regular then T{A\ <g> A 2 ) = T\(A{)T{\ <S> A 2 ). 

(ii) If both A\ and A 2 are strongly regular then T\A\ ® A 2 ) = T\ {A\) 7 2 {A 2 ). 

Remark 2.9. The result of the corollary says that the spectral state factorizes for tensor 
products of strongly regular operators. This corresponds to the formula on page 563 in 
\9$. It should be noticed, however, that this factorisation is in general not valid for tensor 
products of weakly regular operators, since the Laplace transform of a product is not the 
product of the Laplace transforms. We consider below examples of this type. 

Lemma 2.10. Let Oj,6j E Ai, and set {d;a,b} = [d,Tr(a)]*[d,ir(b)]. 

(i) Then one has 

Q(ai <g> a 2 , 6i <8> 62) = T({A; 01,61} ® ir(a* 2 b 2 )) + T(7r(o*6i) <g> {D 2 ; a 2 , b 2 }) . 

(ii) Ifir(a*bi), i = 1,2, are strongly regular then 

Q{ ai ®a 2 M®b 2 ) = T({D 1 ;a 1 ,b 1 }®l)T 2 (n(alb 2 ))+Ti(n(a* 1 b 1 ))T{l®{D 2 -,a 2 ,b 2 }) , 

Proof. Since [D,Tr(ai <g> a 2 )] = [Di,ir(ai)] <g> 7r(a 2 ) + X 7r (ai) ® [D 2 ,ir(a 2 )] we get 

[At(oi®02)]*[At(6i®62)] = [A,7r(ai)]*[A,7r(fei)](2)7r(^6 2 ) 

+tt(o5&i)<8 [D 2 ,7r(a 2 )]*[D 2 ,n(b 2 )\ 
+[A,7r(ai)]*x7r(6i)(8)7r(a2)*[A,7r(6 2 )] 
+X7r(a;)[A,ir(6i)] ® [A,vr(a 2 )]*7r(6 2 ) 

Since Tr (x^4) = for any odd operator A E B($)x) the contributions of the last two lines 
vanish under the trace, hence under T, and we get the first statement. Point (ii) follows 
from Corollary E3J P □ 
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3 The spectral triple associated with a stationary Bratteli 
diagram 

An oriented graph Q = (V, £) is the data of a set of vertices V and a set of edges £ with 

two maps, £ = % V , one assigning to an edge e its source vertex s(e) and the second 

assigning its range r{e). The graph matrix of Q is the matrix A with coefficients 
equal to the number of edges which have source v and range w. 

We construct a spectral triple from the following data (see Figure ?? for an illustration 
of the construction) : 

1. A finite oriented graph Q = (V,£) with a distinguished one-edge-loop I*. We 
suppose that the graph is connected in the stronger sense that for any two vertices 
v\ , v 2 there exists a path from to v\ to V2 and a path from vi to v\ . This is 
equivalent to saying that the graph matrix A is irreducible. We will even require 
that A is primitive (see below). 

Alternatively we can pick a distinguished loop made of p > 1 edges, and resume the 
case described above by replacing the matrix A by A p and considering its associated 
graph Q p instead of Q. 

2. A function f : £ — > £ satisfying: for all e G £ 

(a) if r(e) is the vertex of /* then f (e) = /*, 

(b) if r(e) is not the vertex of I* then f(e) is an edge starting at r(e) and such 
that r(f(e)) is closer to the vertex of /* in Q. 

3. A symmetric subset % = H(Q) of £ x £: 

U C {(£,£') G £ x £ : e ^ e', s(e) = a(e')} 

This can be understood as a graph with vertices £ and edges %, which has no 
loops. We fix an orientation of the edges in %, and write % = % + U %~ for the 
decomposition into positively and negatively oriented edges. 

4. A real number p £ (0, 1). 

r 

Notation. We still denote the range and source maps by r,s on %: % j £ . 
We allow compositions with the source and range maps from £ to V, which we denote 

r 2 ,s 2 

by s 2 ,r 2 ,sr, rs: % V . We denote by H n {Q), or simply by Ii n if the graph is 

understood, the set of paths of length n over Q: sequences of n edges e± ■ ■ ■ e n such that 
r(ffj) = s(ej + i). We also set IIo(^) = V. We extend the range and source maps to paths: 
if 7 = £i---e n then r(j) := r(e n ),s(j) := s(ei). Also, given 7 = e\ ■ ■ ■ £j ■ ■ ■ e n , we 
denote by 7^ = £j the i-th edge along the path. 

The number of paths of length n starting from v and ending in w is then A n vw . We 
require that A is primitive: 3N G N, Vu, u;, > 0. Under this assumption, A has 
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a non-degenerate positive eigenvalue A PF which is strictly larger than the modulus of 
any other eigenvalue. This is the Perron-Frobenius eigenvalue of A. Let us denote by L 
and R the left and right Perron-Frobenius eigenvectors of A (i.e. associated with A PF ), 
normalized so that 

J>, = 1, J>,L, = 1. (8) 

3 3 

Let us also write the minimal polynomial of A as ha(X) = IT? =1 (A — Xk) mk with Ai = 
A PP , m\ = 1. Then from the Jordan decomposition of A one can compute the asymptotics 
of the powers of A as follows: 

i^^A^ + fpfW^, (9) 

k=2 

(a) 

where P^, is a polynomial of degree nik — 1 if n > mk, and of degree less than m& — 1 
if n < nik- 

Let Mj be the algebraic multiplicity of the j-th eigenvalue of A (hence M\ = 1). Let 
RP', for 1 < j < p and 1 < I < Mj, be a basis of (right) eigenvectors of A: AR 3,1 = XjRP' . 
Let also LP' be a basis of left eigenvectors of A normalized so that RP' ■ LP ,l = 5jj>5u>. 
So R 1,1 = R and L ' = L as defined in equation j9]). For 1 < j < p, let us define 

Mj 
J 1=1 v£V 

hen 

Given Q we consider the topological space Hqo of all (one-sided) infinite paths over Q 
with the standard topology. It is compact and metrisable. The set IIoo* of infinite paths 
which eventually become I* forms a dense set. 

Remark 3.1. This construction is equivalent to that of a stationary Bratteli diagram: 
this is an infinite directed graph, with a copy of the vertices V at each level n > 0, and a 
copy of the edges £ linking the vertices at level n to those at level n + 1 for all n (there 
is also a root, and edges linking it down to the vertices at level 0). So for instance, the 
set Tl n of paths of length n, is viewed here as the set of paths from level down to level 
n in the diagram. See Figure [1] for an illustration (the root is represented by the hollow 
circle to the left). 

Given f we obtain an embedding of II^ into n„-(_i by E\ ■ ■ ■ e n i— V e\ ■ ■ ■ s n r(e n ) and 
hence, by iteration, into IIoo* C Hoo- We denote the corresponding inclusion U n <^-> IIoo 
by r. 

Given T~L we define horizontal edges H n between paths of Tl n , namely (7,7') € H n if 
7 and 7' differ only on their last edges e and e', and with (e, e') 6 T~L. For all n, we carry 
the orientation of % over to H n . 
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Figure 1: The stationnary Bratteli diagram associated with the Fibonacci matrix 



The approximation graph G T = (V, E) is given by: 



v = (J v n , v n = r(n n ) c n, 



n>0 



E={jE n , E n = r x r(tf n ) 



n>l 



together with the orientation inherited from H n : so we write E n = E^L)E~ for all n, and 
is 1 = E + U E 1- . Given ft £ H we denote by H n (h), E n (h), and -E^/i) the corresponding 
sets of edges of type h. 

Lemma 3.2. The approximation graph G T = (V,E) is connected if and only if for all 
£,e' G £ with s(e) = s(e') there is a path in TL linking e to e' . Its set of vertices V is 
dense in Hoc- 

Proof. Let x,y £ V, x ^ y, and n the largest integer so that Xi = yi, i = 1, • • • n — 1, 
and x n ^ y n (so s{x n ) = s(y n )). Any path in G T linking x to y must contain a subpath 
linking x n to y n . Hence G T is connected iff the above given condition on % is satisfied. 

The density of V is clear since any base clopen set for the topology of Hoo'. [7] = 
{x € IToo : Xi = rji, i < n}, 7 G n„, n G N, contains a point of V, namely t{^). □ 

Given an edge e G E we write e op for the edge with the opposite orientation: s(e op ) = 
r(e),r(e op ) = s(e). Now our earlier construction [19J yields a spectral triple from the 
data of the approximation graph G T . The C*-algebra is C^oq ), and it is represented on 
the Hilbert space £ 2 (E) by: 



7r(/)V(e) = /«e))V(e) • 



(11) 



The Dirac operator is given by: 



^(e)=p->(e op ) 



e E E. 



'n 



(12) 



The orientation yields a decomposition of i 2 (E) into l 2 (E + ) ®i 2 {E ). 



16 



Theorem 3.3. (C(H OQ (Q)),£ 2 (E), D) is an even spectral triple with Z/2-grading x which 
flips the orientation. Its representation is faithful. Ifri is sufficiently large (i.e. satisfies 
the condition in Lemma \S. #)) then its spectral distance d s is compatible with the topology 
on Hoo(Q)> an d one has: 

d s (x,y) = c xy p n ^ + i b n{x) + b n {y))p n , for x ± y , (13) 

where n xy is the largest integer such that X{ = yi for i < n xy , and b n (z) = 1 if 
T(z n ) ^ z n+ \ and b n (z) = otherwise, for any z £ Iloo. The coefficient c xy is the 
number of edges in a shortest path in "H linking x nxy to y nxy . If rl is maximal: fi = 
{(e, e') £ £ x £ : e ^ e', s(e) = s(e')}, then c xy = 1 for all x,y G IToo. 

Proof. The first statements are clear. In particular the commutator [D,n(f)] is bounded 
if / is a locally constant function. The representation is faithfull by density of V in Iloo. 

If rl satisfies the condition in Lemma 13.21 then the graph G(t) is connected. It is 
also a metric graph, with lengths given by p n for all edges e G E n . By a previous result 
(Lemma 2.5 in |19j ) d s is an extension to Iloo of this graph metric, and as ^ p n < +oo, 
it is continuous and given by equation (|13p (by straightforward generalizations of Lemma 
4.1 and Corollary 4.2 in [IS]). □ 



3.1 Zeta function 

We determine the zeta-function for the triple associated with the above data. 

Theorem 3.4. Suppose that the graph matrix is diagonalizable with eigenvalues Xj, j = 
1, . . .p. The zeta-function £ extends to a meromorphic function on C which is invariant 
under the translation z i— > z + . It is given by 



j=i j 1 



where h is an entire function, and is given in equation (|10p . In particular £ has only 
simple poles which are located at { ~^4e^T~ '■ k £ 1*, j = ■ ■ ■ p} with residui given by 

Res( C , l0SAj + 27Ufc ) = . (14 ) 

- log p - log p 

In particular, the metric dimension is equal to sq = z^p- ■ 
Proof. Clearly 

n>\ 
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The cardinality of E n can be computed by summing, over all vertices v at level n — 1 
and all edges h E T-L with s 2 (h) = v, the number of paths from level down to v at level 
n — 1: 

#E n = Y,Y1 A v7 2 \h) ■ 



v&V h&H 

1 k 



Now since A is diagonalizable, the polynomials Pj in equation ([9]) are all constant and 
can be expressed in terms of its (right and left) eigenvectors R k ' 1 , L k ' 1 ,1 < I < Mj.. 
(These vectors were so normalized that the R k ' 1 are the columns the matrix of change 
of basis which diagonalizes A, and the vectors L k,t are the rows of its inverse.) So from 
equation ([9]) we get for all n, ffE n = 5^fe=i ^k > an< ^ nence 

fc=l n>l k=l kl 

Hence ^ has a simple pole at values z for which /O z Afc = 1, A; = 1, . . .p. The calculation 
of the residui is direct. □ 

The periodicity of the zeta function with purely imaginary period whose length is only 
determined by the dilation factor is a feature which distinguishes our self-similar spectral 
triples from the known triples for manifolds. Note also that £ may have a (simple) pole 
at 0, namely if 1 is an eigenvalue of the graph matrix A. 

Remark 3.5. In the general case, when A is not diagonalizable, it is no longer true that 
the zeta-function has only simple poles. Here the polynomials Pi 3 (n) in equation ([9]) are 
non constant (of degree — 1 > for k = 2, . . .p), and give power terms in the sum 
for C(z) written in the proof of Theorem \3.4\ (one gets sums of the form Xm>i ""A^p™ 2 
for integers a < — 1). In this case ((z) has poles of order m,j at z = logAj lc ^ 2 " fc - 

3.2 Heat kernel 

We derive here the asymptotic behaviour of the trace of the heat kernel Tr (e~ tD2 ) around 
0. We give an explicit formula when the graph matrix A is diagonalizable, and explain 
briefly afterwards (see Remark 13. 7p how the formula has to be corrected when this is not 
the case. 

For r > 0, K(a) > and s£lwe define 

p( W )= £ r(a + ^)e 2 ^. (15) 

k=— oo 

In particular, 

Pr,a(° :=-p r,-,- 16 

' iy ry ry 

is a periodic function of period r, with average over a period p r)(I = ^r(^). 
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Theorem 3.6. Consider the above spectral triple for a Bratteli diagram with graph matrix 
A and parameter p 6 (0,1). We assume that A is diagonalizable with eigenvalues Xj, 
j = 1, . . .p. Then the trace of the heat kernel has the following expansion as t — » + ; 



Tr(e- tD )= £ C H P-2iogp.log A, 

j:\\j\>l 



-log* 
-21ogp 



+ h(t) (17) 



where is given in equation (|10p . jo is the index for the eigenvalue Xj = 1 (setting 
=0 if no eigenvalue equal to 1 is present) and h is a smooth function around 0. In 
particular, the trace of the heat kernel is equivalent to 

j-r-,2 t— ifl i log App 

Tr(e- tD ) ~ Cip_2iogp,iogA PF (-logt) t^T . 

Proof. From equation Q and the proof of Theorem 13.41 we have, for all n, ifE n = 
£^ =1 C^A«,andso 



£#£, 



-t P - 



We split the sum into several parts: the one over I = {j : \Xj\ < 1}, over J = {j : \Xj\ > 
1}, and possibly the terms for the eigenvalues of modulus 1. Clearly 



££4^" 

jei n 



<£i4i£N n < 

jel n 



3G 



so the eigenvalues of modulus smaller than 1 do not contribute to the singular behaviour 
at t = 0. 

To treat the eigenvalues of modulus larger than 1 (with indexes in J), set v = p~ 2 
and ctj = _2k)gp (principal logarithm). Note that 9?(ay) > for j € J. Then 



jGJ n 



£ °« £ ^ «p(-tw») 
j'eJ n=0 



£4 

J6J 



+oo 



£ w nc ^ exp(-to") - £ w"" * exp(-to-") 



n=l 



The term ""j exp(— tv n ) being bounded at we only need to concentrate on 

the first sum. Clearly V*4 En^-oo ^ exp(-iu n ) = /(p§A) where 



> log V ' 



n+s\ 
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By standard arguments the series is uniformly convergent defining a smooth 1-periodic 
function /. It follows that the singular behaviour of ^^oV nQ!j ' exp(— as t — > + is 
given by f(logt)t~ a i . Since the Fourier coefficients of / are given by 



oo 



r(a,- + 



2mk i 

e los ^ x exp(-v x )e 2mkx dx = - — - logv ' 

oo log" 



we see that f(s) = _ 21 1 og/3 p(-2 log p, a>j,-s). 

Consider now an eigenvalue of modulus 1, so ay = up for some tp € (— ir, ir]. We claim 
that sup t | Y,n=o v m ^e~ tvn | < oo if and only <p ^ 0. Notice that sup t | £n=o v m ^e~ tvn \ < 
oo if and only if sup i6N iWe [o,i) I Yln=-L v mv e~ vn+ " \ < oo. We split the sum 



-n-t-w j 

n=0 n=0 n=0 



and replace the second sum by an integral (x — > 1 — e v " " is strictly decreasing and 
finite over R): 

l i r w i-e-y , i r 1 - e -» 







(1 - e v )dx = / dy < / dy < oo. 

logv J v -l+u, y log Wo V 



Now sup^ | X^n=o v~ mLp | is finite whenever <p ^ and so the eigenvalues of modulus 1 
which are not equal to 1 do not contribute to the singular behaviour at t = 0, either. 

Finally suppose that one of the eigenvalues, say A, , is equal to 1. Then its contribu- 
tion to the trace of the heat kernel is ^ n e~ tp 2 ™ which we may study through the 
integral: 

e-f-ck = J- f° e —ds ^° " l0g * - bgt 



log v J t s log v 2 log p 

Putting all contributions together we obtain the statement. □ 

We could have determined asymptotic expansion of the trace of the heat kernel using 
the inverse Mellin transform of the function Q(2s)T(s). This is, of course, a lot more 
complicated than the direct computations. When using the inverse Mellin transform of 
£(2s)r(s) one can see that the origin of the periodic function p is directly related to the 
periodicity of the zeta function and that the appearence of the term log t in the trace of 
the heat kernel expansion arises from a simple pole of at s = which amounts to a 
double pole of £(2s)r(s) at s = 0. 

Remark 3.7. In the general case, when A is not diagonalizable, the formula in The- 
orem 13.61 has to be modified. This correction affects however only those terms which 
correspond to eigenvalues whose algebraic multiplicity is different from the geometric 
multiplicity. In particular, the asymptotic behaviour as t — > + remains unaffected. 
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We just indicate how things are altered. Let rrij be the size of the largest Jordan block 
of A to eigenvalue Xj and Pj the polymomial of degree rrij — 1 as above in equation ([9]) . 
The contribution of eigenvalues of modulus greater than one in equation (|17p now reads 

j:\X d \>l ' 

where Sj = °f g p ■ While the contribution from Xj Q = 1 becomes now \C^Pj (^^) , 

where Pj is a "primitive" of Pj , namely if Pj = ^2 akX k then Pj = ^ 2 ^ +] ^ 1 afcX fc+1 , 
But unlike the diagonalizable case, there is also a singular contribution from the other 
eigenvalues Xj of modulus 1 not equal to 1: \Xj\ = l,Xj ^ 1. 

In view of Remark \3.5\ we also notice that the poles of order rrij in the zeta function 
at z = Sj + ^\og p 5* ue r i se t° the log-terms Pj (j^p-), i-£- to terms log(i) 9 for q < rrij — 1 
(the degree of Pj). In general a pole of order m + 1 in the zeta function gives rise to 
log-terms log(t) q , q < m. 

3.3 Spectral state 

There is a natural Borel probability measure on IIoo. Indeed, due to the primitivity of 
the graph matrix there is a unique Borel probability measure which is invariant under 
the action of the groupoid given by tail equivalence. We explain that. 

If we denote by [7] the cylinder set of infinite paths beginning with 7, then invariance 
under the above mentionned groupoid means that ^([7]) depends only on the length I7I 
of 7 and its range, i.e. ^([7]) = /u(|7|, r (l))- By additivity we have 

= £ MM 

e:s(e)=r(-y) 

which translates into 

fi(n, v) = A vw [i{n + l,w). 

w 

The unique solution to that equation is 

fi(n, v) = X~™R V 

where R is the right Perron- Frobenius eigenvector of the adjacency matrix A, normalized 
as in equation ([8|). So if 7 E II n is a path of length n, then ^([7]) = Ap"i? r ( 7 ). 

Theorem 3.8. All operators of the form ir(f), f £ C^IIoo), are strongly regular. More- 
over, the measure fi defined on f G C(Jloo) by 

M/) := rw)) 

is the unique measure which is invariant under the groupoid of tail equivalence. 
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Proof. Let / be a measurable function on Iloo, and set 

ff ,_ ^EMf)\E n ) _ £ egJ E„ /(g(<0) 

To check that the sequence {p n {f))n converges it suffices to consider / to be a char- 
acteristic function of a base clopen set for the (Cantor) topology of Iloo. Let 7 be a 
finite path of length I7I < n and denote by x-y the characteristic function on [7]. Then 
X 7 (s(e)) is non-zero if the path s(e) starts with 7. Given that the tail of the path s(e) is 
determined by the choice function r, the number of e £ E n for which % 7 (s(e)) is non-zero 
coincides with the number of paths of length n — \ — 1 which start at r{pi) and end at 
s(s(e) n -i) = s 2 (h). Hence 



Ex7(-(e))=E^r l7l_l 

eeEn hen 



r (l) s 2 (h) 



As noted before in the proof of Theorem 13.41 the cardinality of E n is asymptotically 
C-^App, so we have 

Mn(X 7 ) = = A pp -Q- R r{l) L s*{h) (1 + o(l)) • 

ft n n hen 

Set U v = (\pp /Cu) YliheH RvL s 2( h -)). We readily check that U is a (right) eigenvector of 
A with eigenvalue A PF , and since its coordinates add up to 1, we have U = R. So we get 

A*n(x 7 ) = Api 7l i? r(7) (l + o(l)) ?W+00 > Ap] 7 '^,,^) = ^([7]) . 
Now Corollary 12.51 implies that vr(/) is strongly regular and 7~(7r(/)) = □ 

Below, when we discuss in particular the transversal part of the Dirichlet forms for a 
substitution tiling, we will encounter operators which are weakly regular, but not strongly 
regular, and so, a priori, the spectral state is not multiplicative on tensor products of 
these. The following lemma will be useful in this case. 

Lemma 3.9. Consider the above spectral triple for a Bratteli diagram with graph matrix 
A and parameter p £ (0, 1). Let A £ B(Sy) be such that A n n ~°° e mip for some if £ (0, 2ir). 
Then 

, n 2 /— - £ Iog - log App 

Tr(e-* D ^) ~ e 4 2i° g pC w p_ 21ogPilogA p P+ ^(-logt)t 2io gP . 
In particular, A is weakly regular and T(A) = 0. 

Proof. As before we set a = j^;, v = p~ 2 , a = !?f logp so ^ na ^ we Fiave to de- 
termine the asymptotic behaviour of Ylri=i e m{p v a ^ n+a ^ e~ v " +rT when a — > —00. Since 
^2n=-oo e m{p v a ( n+a ') e~ vn+ ° is absolutely convergent the sum over N has the same asymp- 
totic behaviour than the sum over Z. Now 

00 00 



E e mip v a ( n+cr ^e~ vn+ ' T = e~ lcrip ^ 



- log A PF + up 
pQogv, - , -a) 



log v log V 
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from which the first statement follows. T(A) is equal to the mean of the function a i— > 
Qoet ?; logApF+ 



e mLp j-i— p(log v, log \oIv~ i(f ' ~~ a ) w hi cn i s zero since for all integer k is tp + 27r£; 7^ 0. □ 



We now consider the tensor product of two spectral triples associated with Bratteli 
diagrams, one with parameter p < 1 and Perron Frobenius eigenvalue A PP the other with 
parameter p' < 1 and Perron Frobenius eigenvalue A PF . 

Lemma 3.10. Let A be as in the last lemma, then T{A ® 1) = if for all integer k, k' 

09 + 2irk + 2-K-^rk' / 0. (18) 
log p' 

Proof. By the above general results T(^4 ® 1) is equal to the mean of the function s 1— > 
J&pQogv, log f P g F +^ , ^)p(log^, ^) devided by the mean of the function 

s t-> p(log v, nllf 1 , I ^)p(logw / , 1^7). The latter is always strictly greater than 

zero. By developing the two p-functions into Fourier series one sees that the former mean 

^ sip 2nks „ 2-Kk 1 s 

is zero if s 1— > e lo s K e lo E"e lo s"' is oscillating for all k, k . □ 

Motivated by this lemma we define: 

Definition 3.11. Let p,p' € (0,1). We call p € (0,2-/r) a non resonant phase for (p, p') 
if 

ip + 2-Kk + 2-K-^-k' ^ 0, Vk,k'eZ. (19) 
log // 

If is irrational and for instance c/? = 2it ^4 , then 99 is resonant and, as follows 
easily from the calculation in the proof, T(A ® 1) 7^ 0. In particular, T(^4 (g> 1) 7^ 
Ti(A)T 2 (l). 

3.4 Quadratic form 

In Section 12.31 we considered a quadratic form which we now specify to the context of 
spectral triples defined by Bratteli diagrams. Let £ and sq be the zeta function with its 
abscisse of convergence and p be the spectral measure on Iloo as defined in Section 13.11 
and 13.31 Equations (fTT|) and (|12p imply, for e £ E n 

[D,7r(f)Me) = (S«f)p(e°v), with 5 e f = /(r(e)) ~ f{s{e)) . (20) 



P" 



Equation ([6]) therefore becomes 



Q(f,g) = ]ha-^Tr(\Dr{DMf))*{DM9)]) 



lim + 7rrE#^' ons in(f,g) (21) 
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with 

q n (f,g) = -± r £ 6J5 e g. (22) 

* En ett 

Notice that lim a _ +a + Xm>i i^^nP ns = 1< We thus have the following simple result: 

Lemma 3.12. Q(f,g) = lim^-nx, q n (f, g) provided the limit exists. 

The following can be said in general. 

Proposition 3.13. The quadratic form Q is symmetric, positive definite, and Markovian 
on V = {f e i£Cnoo, dp) : Q(fJ)< +00} . 

Proof. The form is clearly symmetric and positive definite. Let e > and consider a 
smooth approximation g e of the restriction to [0, 1]: g e {t) = t for t £ [0, 1], — e < g e (t) < 
1 + e for t E 1, and < g t (t) — <7e(0 <t — t' for all t < t' . This last property implies that 
l^e 9e°f \ < l^e/| so that Q(g e ° f, g e ° f) < Q(fyf)- This proves that Q is Markovian. □ 

For a precise evaluation of this quadratic form on certain domains we have, however, 
to consider more specific systems. 

3.5 A simple example: The graph with one vertex and two edges 

Remember that one difficulty with the Dirichlet form defined by our spectral triple (equa- 
tion ([6]) in Section I2.3() is that we need to specify a core for the form. We provide here 
an example where such a core can be suggested with the help of an additional structure. 

We consider the graph Q which has one vertex and two edges, which are necessarily 
loops. Call one edge and the other 1. Let I* be the edge 0. Then it is clear that 
IIoo can be identified with the set of all {0, l}-sequences and IIoo* with those sequences 
which eventually become 0. We consider the spectral triple of Section [3] associated with 
the graph Q and parameter p = ^. There is not much choice for the horizontal egdes: 
~H = {(0, 1), (1,0)}, nor for f : r(l) = and f(0) = 0. We will look at this system from 
two different angles justifying two different cores for the Dirichlet form. 

Group structure The space of {0, l}-sequences carries an Abelian group structure. 
In fact, if we identify n n with using 7 1— > Yli=i 7«2 l_1 then we can write IIoo as the 
inverse limit group 



IIoo = lim I Z 2 n ^— 



where 7r™ +1 maps m mod 2 n+1 to m mod 2 n . 

It follows that C(Hoo) is isomorphic to the group C*-algebra of the Pontrayagin dual 



IIoo = lim 



1/21+1 



This suggest that a reasonable choice for the domain of the form would be the (algebraic) 
group algebra CIIoo whose elements are finite linear combinations of elements of IIoo ■ Now 
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CIloo corresponds to the subalgebra Ci c X^oo) C C(Hoo) of locally constant functions. 
We have lim n _ 5 . 00 q n (f, g) = if / and g are locally constant|f| and so, by Lemma [3.121 Q 
exists on that core. However, Q is identically on the domain defined by the core. So 
the algebraic choice of the domain which was motivated by the group structure is not 
very interesting. 

Embedding into S 1 Note that IL n can be identified with via the map 7 1— > 

exp(2™^ =l7i 2-J +1 (-)) so that LL^ = ftoo = {exp(2mp(-)) : p € Z[±] n [0,1)}. This 
suggest to view Hcq* (via 7 1— > exp(2-7^^ X^i>i 7-;2~* +1 ) G 5" 1 ) as a dense subset of 5 1 . 
The inclusion LToo* 5* 1 extends to a continuous surjection 77 : — > S 1 which is 
almost everywhere one-to-one. Furthermore the push forward of the measure p on Iloo 
is the Lebesgue measure on S 1 . Hence the pull back 77* induces an isometry between 
L 2 ^ 1 ) and L 2 (n M ). It follows that ^(C^S 1 )) is dense in L 2 (noo). This suggests to 
take as core for the quadratic form the pull backs of trigonometric polynomials over 
S l . Now one sees from equation (I22p that lim n _ 5 . 00 q n (f, g) exists, provided p = 5. In 
fact, linin^oo q n {f,g) = §f ), showing that Q has infinitesimal generator equal to the 
standard Laplacian on S . 



3.6 Telescoping 

There is a standard equivalence relation among Bratteli diagrams which is generated 
by isomorphisms and so-called telescoping. Since we are looking at stationary diagrams 
we consider stationary telescopings only. Then the following operations generate the 
equivalence relation we consider: 

1. Telescoping: Given the above data built from a graph Q = (V, £), and a positive 
integer p, we consider a new graph Q p := (V P ,S P ) with the same vertices: V p = V, 
and the paths of length p as edges: £ p = H p (Q). The corresponding parameter is 
taken to be p p = pP . 

2. Isomorphism: Given two graphs as above Q = (V, £), Q' = (V',£'), we say that 
the corresponding stationary Bratteli diagrams are isomorphic if there are two 
bijections V — > V, £ — > £' which intertwine the range and source maps. We need 
in this case that the associated parameters be equal, and the sets of horizontal 
edges isomorphic (through a map which intertwines the range and source maps). 

We show now that this equivalence relation leaves the properties of the associated spectral 
triple invariant: 

(i) The zeta functions are equivalent, so have the same spectral dimension: so = Tr^p- ; 

(ii) The spectral measures are both equal to the invariant probability measure p on 



3 If / is locally constant, then S e f = in equation (1201) . for all e G E n for n large enough. 
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(iii) Both spectral distances generate the topology of 11^ (provided H is large enough 
as in Lemma l3,2p . and are furthermore Lipschitz equivalent. 



The invariance under isomorphism is trivial. We explain briefly how things work under 
telescoping. The horizontal edges for Q p are given as for Q, by the corresponding subset 

W C {(£,£') G £ p x £P : e £ e', s(e) = s{e')} , 

namely: writting e = £\ . . . e p and e' = e' x . . . s 1 (by definition £ v = U p (Q)), there is an 
1 < i < p such that £j = e'- for 1 < j < i, and (e^, e'j) E So we have the identifications 

p—i p— l 

H p = [^J H np +i , E p = E n p+i 

i=0 i=0 

which allows us to determine the approximation graph G? = (V P ,E P ), and yields a 
unitary equivalence i 2 (E) = £ 2 (E P ). We identify the two Hilbert spaces i 2 (E) = £ 2 (E P ) 
and the representations ir = ir p , while the Dirac operators satisfy: 



, [ f(E) -> l 2 (E p ) 

Dp = W> D W with V, .-WU.^.wW.n-* m od P . 



From the inequalities 1 < < we deduce that the zeta functions are equivalent, 

and that both spectral triple have the same spectral dimension so, and give rise to the 
same spectral measure /i. By Theorem 13.31 both Connes distances generate the topology 
of IIoo = n^o, provided % is large enough. Let us call dP s the spectral metric associated 
with G p , with corresponding coefficients n^c^,^ as in equation (|13p . Writing n 
pn xy + k xy , for some k xy < p — 1, we have 



XII 



p-1 

d s (x,y) = c xyP k ^{p p ) nl y + (p p ) n Y,( bnp + k ^ + b np +k(y))p k , 

n>n p xy k=0 

Now we see that b np+ k{z) = 1 =>- b P l (z) = 1, while if b np+ k{z) = for all k = 0, ■ • -p — 1, 
then 6^(z) = too, so that one has fe^(z) < Y^=o b np +k( z ) < pbn(z). We substitute this 
back into the previous equation to get the Lipschitz equivalence: 

c p d p (x,y) <d s {x,y)<pp p C p d p (x,y), 

with c p , C p , the respective min and max of c xy j(? xy (which only depends on 7i and p). 



4 Substitution tiling spaces 

Bratteli diagrams occur naturally in the description of substitution tilings. The path 
space of the Bratteli diagram defined by the substitution graph has been used to describe 
the transversal of such a tiling |11U17| . As we will first show, an extended version can also 
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be used to describe a dense set of the continuous hull f2 CT of the tiling and therefore we will 
employ it and the construction of the previous section to construct a spectral triple for 
Q a . We then will have a closer look at the Dirichlet form defined by the spectral triple. 
Under the assumption that the dynamical spectrum of the tiling is purely discrete we 
can identify a core for the Dirichlet. We can then also compute explicitely the associated 
Laplacian. 



4.1 Preliminaries 

We recall the basic notions of tiling theory, namely tiles, patches, tilings of the Euclidean 
space M. d , and substitutions. For a more detailed presentation in particular of substitution 
tilings we refer the reader to |14| . A tile is a compact subset of M d that is homeomorphic 
to a ball. It possibly carries a decoration (for instance its collar). A tiling of M. d is a 
countable set of tiles (ti)«e/ whose union covers M. d and with pairwise disjoint interiors. 
Given a tiling T, we call a patch of T, any set of tiles in T which covers a bounded and 
simply connected set. A prototile (respectively protopatch) is an equivalence class of tiles 
(respectively patches) modulo translations. We will only consider tilings with finitely 
many prototiles and for which there are only finitly many protopatches containing two 
tiles (such tilings have Finite Local Complexity or FLC). 

The tilings we are interested in are constructed from a (finite) prototile set A and a 
substitution rule on the prototiles. A substitution rule is a decomposition rule followed 
by a rescaling, i.e. each prototile is decomposed into smaller tiles, which, when stretched 
by a common factor 9 > 1 are congruent to some prototiles. We call 9 the dilation factor 
of the substitution. The decomposition rule can be applied to patches and whole tilings, 
by simply decomposing each tile, and so can be the substitution rule when the result of 
the decomposition is stretched by a factor of 9. We denote the decomposition rule by 5 
and the substitution rule by a. In particular we have, for a tile t, 6(t + a) = 5(t) + a and 
o~(t + a) = a(t) + 6a for all a € M. d . See Figure [2] for an example in R 2 . 



Figure 2: A substitution (chair tiling): decomposition and dilation. 



A patch of the form o~ n (t), for some tile t, is called an n-supertile, or n-th order 
supertile. A rescaled tile 6 n t will be called a level n tile but also, if n = —m < 0, an 
m-microtile, or m-th order microtile. 
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A substitution defines a tiling space Q a : the set of all tilings T with the property 
that any patch of T occurs in a supertile of sufficiently high order. 

We will assume that the substitution is primitive and aperiodic: there exists an integer 
n such that any n-supertile contains tiles of each type and all tilings of Q a are aperiodic. 
This implies that by inspection of a large enough but finite patch around them the tiles 
of £l a can be grouped into supertiles (one says that a is recognizable) so that 5 and a 
are invertible. In particular, a is a homeomorphism of CT if the latter is equipped with 
the standard tiling metric [I]. 

We may suppose that the substitution forces the border |17j . The condition says that 
given any tile t, its n-th substitute does not only determine the n-supertile o~ n (t), but 
also all tiles that can be adjacent to it. This condition can be realized for instance by 
considering decorations of each types of tiles, and replacing A by the larger set of collared 
prototiles. 

There is a canonical action of M. d on the tiling space Q a , by translation, which makes 
it a topological dynamical system. Under the above assumptions, the dynamical system 
(J7 CT ,R d ) is minimal and uniquely ergodic. The unique invariant and ergodic probability 
measure on Q a will be denoted u. 

A particularity of tiling dynamical system is that they admit particular transversals to 
the Reaction. To define such a transversal S, we associate to each prottoile a particular 
point, called its puncture. Each level n tile being similar to a unique proto-tile we may 
then associate to the level n tile the puncture which is the image of the puncture of 
the proto-tile under the similarity. The transversa]^ S is the subset of tilings T £ ft 
which has the puncture of one of its tiles at the origin of The measure u induces an 
invariant probability measure on H which gives the frequencies of the tiles and patches. 

4.2 Substitution graph and the Robinson map 

The substitution matrix of the substitution a is the matrix with coefficients equal to 
the number of tiles of type t{ in o~(tj). The graph Q of Section [3] underlying our construc- 
tions will be here the substitution graph which is the graph whose graph matrix is the 
substitution matrix. More precisely, its vertices v € V are in one-to-one correspondence 
with the prototiles and we denote with t v the prototile corresponding to v € V. Between 
the vertices u and v there are A uv edges corresponding to the A uv different occurrences 
of tiles of type t u in a(t v ), see Figure [3j Here we call u (or t u ) the source, and v (or t v ) 
the range of these edges. Notice that the Perron-Frobenius eigenvalue of A is the d-ih 
power of the dilation factor 6: A PF = 9 d . The asymptotics of the powers of A are given 
by equation Q as before. The coordinates of the left and right Perron-Frobenius eigen- 
vectors L, R, normalized as in (JBJ) are now related to the volumes and the frequencies of 
the prototiles as follows: for all v € V we have 

freq(tu) = R v , vol(t„) = L v (23) 
4 Sometimes S is referred to as the canonical transversal or the discrete hull. 
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Figure 3: An illustration of the Bratteli diagram associated with the chair substitution. 



where freq(i„) is the frequency and vol(t v ) the volume of t v , the volume being normalized 
as in equation (JS]) so that the average volume of a tile is 1. 

Given a choice of punctures to define the transversal H of Q a there is a map 

onto the set of half-infinite paths in Q. Indeed, given a tiling T 6 3 (so with a puncture 
at the origin) and an integer n E N, we define: 

• v n (T) € V to be the vertex corresponding to the prototile type of the tile in a~ n (T) 
which contains the origin; 

• e n (T) E £ to be the edge corresponding to the occurrence of v n -\{T) in a(v n (T)). 

Then 1Z(T) is the sequence (e n (T)) n >i. We call 7Z the Robinson map as it was first 
defined for the Penrose tilings by Robinson, see |14j . 

Theorem 4.1 (|17|). 1Z is a homeomorphism. 

We extend the above map 1Z to the continuous hull il a . The idea is simple: the 
definition of TZ makes sense provided the origin lies in a single tile but becomes ambiguous 
as soon as it lies in the common boundary of several tiles. We will therefore always assign 
that boundary to a unique tile in the following way. 
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We suppose that the boundary of the tiles is sufficiently regular such that there exist 
a vector v G M. d such that for all points x of a tile t either 3i] > 0, Ve G (0, rf) : x + ev G t 
or 3r] > 0, Ve G (0,7/) : x + ev ^ t. This is clearly the case for polyhedral tilings. We fix 
such a vector v. Given a prototile t (a closed set) we define the half-open prototile [t) as 
follows: 

[t):= {x Et : 3t] > Ve G [0, rf) : x + ev G t} . 

It follows that any tiling T gives rise to a partition of R rf by half-open tiles. We extend 
the Robinson map to 

where rL^+oo is the space of bi-infinite sequences over Q using half-open proto-tiles as 
follows. For neZwe define 

• v n (T) G V to be the vertex corresponding to the prototile type of the half open tile 
in a~ n (T) which contains the origin. So v n (T) corresponds to 

— the n-th order (half-open) supertile in T containing the origin, for n > 0, 

— n-th order (half-open) microtile in 5~ n (T) containing the origin, for n < 0; 

• £ n (T) G £ to be the edge corresponding to the occurrence of v n -\{T) in a(v n (T)). 

And we set 1Z(T) to be the bi-infinite sequence 1Z(T) = ( £ n(T')) ngZ - 

Remark 4.2. As in Remark 13. II we can see this construction as a Bratteli diagram, but 
the diagram is bi-infinite this time. There is a copy of V at each level n G Z and edges 
of £ between levels n and n + 1. Level corresponds to prototiles, level 1 to supertiles 
and level n > 1 to n-th order supertiles, while level —1 corresponds to microtiles and 
level n < — 1 to n-th order microtiles. For the "negative" part of the diagram, we can 
alternatively consider the reversed substitution graph Q = (V, £) which is Q with all 
orientations of the edges reversed. The graph matrix of Q, is then then transpose of the 
substitution matrix: A = A T . So for n < 0, there are edges of £ between levels n and 
n — 1: there are A uv = A vu such edges linking u to v. 

As for Theorem 14. H one proves, using the border forcing condition, that 1Z is injective. 

Given a path £ G IL-^+oq, and m < n G ZU {±oo}, we denote by £[ m , n ]£(m,n]£[m,n) 
and i(jn )n ), its restrictions from level m to n (with end points included or not). Also £ n 
will denote its n-th edge, from level n to level n + 1, n 6 Z. We similarly define Ti m) n 
(with end points included). For instance IIo^oo is simply = Il 00 (^). 

We say that an edge e G £ is inner if it encodes the position of a tile t in the supertile 
p such that 3n > 0,Ve G [0,r/) : t + ev G p. This says that the occurence of t in p does 
not intersect the open part of the border of p. 

It is not true that 7Z is bijective but we have the following. 

Lemma 4.3. X := \mlZ contains the set of paths £ G IL-o^+oo such that £(-0^0] has 
infinitely many inner edges. 
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Proof. Recall the following: If (A n ) n is a sequence of subsets of M. d such that A n+ i C A n 
and diam(^4 n ) — > then there exists a unique point such that x € f] n A n . 

By construction, for any tiling T, one has € (~) n <o[v n (T)). Hence ^ = 1Z{T) 
whenever rin<o[' s (^™)) ^ ®> wnere [s(£o)) is the standard representative for the half-open 
prototile of type s(£o) an d [ s (£n)) the half-open n-th order microtile of type s(£ n ) in 
I s (Co)) which is encoded by the path 6 nj ol- 

Suppose that £ n is inner, then [s(£ n )) does not lie at the open border of [r(£ n ) = 
s(£ n+ l)). Hence [s(£ n )) n [s(£ n+ i)) = [s(f n )] n [s(£ n+ i)) where [s(£ n )] is the closure of 
[s(£ n )). Suppose that infinitely many edges of £(_oo,ol are inner, then 

n Kwicfiku)) 

n<0:£„ inner «<0 

showing that the r.h.s. contains an element, and hence £ € im7£. □ 

Corollary 4.4. T/ie sei X Z5 a dense and shift invariant subset ofYl—^ -^^. 

Proof Shift invariance is clear. Denseness follows immediately from Lemma [4,31 □ 

In particular, for n € N, each element of H- n ,n can be the middle part of a sequence 
in lZ(p, a ), that is, for all 7 £ H- n ,n there exists T E Q a such that 7£(T)r njn i = 7. 

Remark 4.5. For v £ V, let n^.^ ^ be the set of bi-infinite paths which pass through v 
at level 0, and set X v = X Pi ^ Then 7?. yields a bijection between E^ x [t v ) and 
X 1 ', where ^ is the prototile corresponding to v and E^ its acceptance domain (the set 
of all tilings in H which have t v at the origin) . 

Notice that H-ooq can be identified with Hoo(G), where Q is the graph obtained from 
Q by reversing the orientation of its edges: one simply reads paths backwards, so follows 
the edges along their opposite orientations. One sees then, that the Robinson map yields 
a homeomorphism between E^ and +QO = H'^, and a map with dense image from [t v ) 

into n-_ OOi0 = n^(g). 

4.3 The transversal triple for a substitution tiling 

Our aim here is to construct a spectral triple for the transversal E. We apply the general 
construction of Section [3] to the substitution graph Q = (V, £). We may supposqij that 
the substitution has a fixed point T* such that RM d is covered by the union over n of 
the n-th order supertiles of T* containing the origin. It follows that 1Z(T*) is a constant 
path in IL_ 00t+00 (Q), that is, the infinite repetition of a loop edge which we choose to be 
e*. We then fix r, take p = pt r as a parameter, and choose a subset 

U tr C U(Q) = {(£,£') e £ x £ : e + e', s(e) = s(e')} 

5 This can always be achieved by going over to a power of the substitution. 
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which we suppose to satisfy the conditions of Lemma 13.21 if s(e) = s(e') there is a path 
of edges in % linking e with e' . The horizontal edges of level n G N are then given by 

H tr ,n ={(r/e,r ? e') : r, G U n ^(g), (e, e') G H tr ) C U n (G) x U n (G) . 

They define the transverse approximation graph Gt r , T = (Vt r ,Et r ) as in Section [3] 

V tr = (J Vtr,n , ^r,n = r(II n (g)) C 11^(0) , 

n 

Etr = l^J -Etr,n > ^tr.n = r x T {Htr,n) , 

n 

together with the orientation inherited from TL: so Et r>n = E^, n U E^ rn for all n G N, 
and 2t^ r = 22^ U E^ r . We also write E n (h) = r x T(Ht r , n (h)) where, if /i = (e,e'), then 
Htr,n(h) = {(rje,rj£') : r] G Tl n -i(Q)}. By our assumption on % the approximation graph 
Gtr,r is connected, and its vertices are dense in n oo (t/). 

An edge h G "Htr has the following interpretation: The two paths r(s(/i)) and r(r(/i)) 
have the same source vertex, say Vo, they differ on their first edge and then, at some 
minimal n/j > 0, they come back together coinciding for all further edges. This is a 
consequence of the property of r. Let us denote the vertex at which the two edges come 
back together with Vh- Neglecting the part after that vertex we obtain a pair (7,7') of 
paths of length which both start at vq and end at Vh- Reading the definition of the 
Robinson map 1Z backwards we see that the pair (7,7') describes a pair of tiles (t,t') of 
type vo in an n/j-th order supertile of type Vh- Of importance below will be the vector 
r/j G M, d of translation from t to t' . 

The interpretation of an edge e G Et r ^ n {h) (so an edge of type h) is similar, except 
that the paths r(s(e)) and r(r(e)) coincide up to level n and meet again at level n + n/j. 
In particular e describes a pair of n-th order supertiles (i, t') of type vq in an (n + n/J-th 
order supertile of type Vh- If one denotes by r e G M. d the translation vector between t 
and t' then, due to the selfsimilarity, one has: 

r e = n r h . (24) 

See Figure 3] for an illustration. 

Theorem 13.31 provides us with a spectral triple for the algebra C(LI 00 (t/)). We adapt 
this slightly to get a spectral triple for C(H). Since the ra-th order supertiles of T* on 
eventually cover M. d , 1Z identifies I1^ (C7) with the translates of T* which belong to 5. 
We may thus consider the spectral triple (C(H),.f)£ r , Dt r ) (which depends on pt r and the 
choices for r and %)) with representation and Dirac operator defined as in equations (|lip 
and (Q2]) by: 

Sitr = P{Etr) , M/M<0 = fipr\a{e)))<p{e) , D tr ip(e) = / t 7V(e° p ) , e G E tr>n . 

We call it the transverse spectral triple of the substitution tiling. By Theorem 13.31 it is an 
even spectral triple with grading \ (which flips the orientation). Also, since Titr satisfies 
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Figure 4: A doubly pointed pattern associated with a horizontal arrow e 6 E tr ^(h). The arrow 
represents the vector r e . Here n = 2 (the paths have lengths 2), and rih = 1 (the paths join 
further down at level n + = 3). 

the hypothesis of Lemma 13,21 as noted above, the Connes distance induces the topology 
of E. By Theorems 13,41 and 13,81 the transversal spectral triple has metric dimension 
and its spectral measure is the unique ergodic measure on H which is 
invariant under the tiling groupoid action. 

For v £ V, we will also consider the spectral triple (C (Bt v ) , Sj^ r , Df r ) for = 
IZ^ 1 (ll^ (^)) : the acceptance domain of t v (see Remark 14. 5p . We call it the transverse 
spectral triple for the prototile t v . It is obtained from the transverse spectral triple by 
restriction to the Hilbert space S)'1 r = £ 2 (E% r ) where Ef r are the horizontal edges between 
paths which start on v. This restriction has the effect that 



that is, up to a perturbation which is regular at St r , the new zeta function is Ry = freq(t„) 
times the old one. It hence has the same abscissa of convergene, s\ r = St r , but its residue 
at str is freq(£y) times the old one. 

Like for E the Connes distance induces the topology of E^. Finally the spectral 
measure \j^ r is the restriction to E^ of the invariant measure on S, normalized so that 
the total measure of is 1. The spectral triple for E is actually the direct sum over 
v € V of the spectral triples for E^ . 

We are particularity interested in the quadratic form defined formally by 



Cv = RvC + reg. 



Qtr(f,g) = Ttr([D tr ,TT t r{f)}*[D tr ,7r tr (g)}) 



= lim 

Sistr 



Tr^ r (\D tr \- s [D tr ,7T tr (f)}* [Dtr,7T t r(g)}) 
CtrO) 



(25) 
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We emphasize that this expression has yet little meaning, as we have not yet specified a 
domain for this form. For example, while strongly pattern equivariant functions |18| are 
dense they do not form an interesting core, as Qt r vanishes on such functions (see the 
paragraph on the transversal form in Section 14. 7p . 

4.4 The longitudinal triple for a substitution tiling 

We now aim at constructing what we call the longitudinal spectral triple for the substi- 
tution tiling which is based on the reversed substitution graph Q = (V, £) (Q with all 
orientations of the edges reversed, so with adjacency matrix A = A T ). Set e* = e* and 
choose t. We take p = pi g as a parameter, and choose a subset 

Hl g CH(§) = {{£,?) G £x £ : e^r, s(e) = s(?)} 
= {{e,e')e£x£ : e^e', r{e) = r{e')} 

again satisfying the condition of Lemma 13.21 We denote the horizontal edges of level 
n e N by 

Hi g , n ={(ve,7]?) : V G ^n-i{G), (e,?) € Hi 9 ) C U n (g) x Ti n {Q) , 
and define the longitudinal approximation graph G^ jT = (Vi g ,Ei g ) as in Section [3] by 

Vig = \jvi a , n , v lg , n = T(u n (g))cu* 00 (G), 

n 

Big = U ^ l 9' n ' Eig, n = r X T(Hi gtn ) , 

n 

together with the orientation inherited from 7ii g : so Ei g>n = E^ n U E^ n for all n € N, 
a ndE lg = E+UE lg . 

With these choices made, Theorem 13.31 provides us with a spectral triple for the 
algebra C(n oo (a)). 

A longitudinal horizontal edge h G Hi g has the following interpretation: As for the 
transversal horizontal edges, r(s(/i)) and r(r(h)) start on a common vertex vq, differ on 
their first edge and then come back to finish equally. To obtain their interpretation it is 
more useful, however, to reverse their orientation as this is the way the Robinson map 1Z 
was defined. Then h = (s,£ f ) with r(e) = r(e') determines a pair of microtiles (t, t') of 
type s(e) and s(e'), respectively, in a tile of type r(e). The remaining part of the double 
path (t(e),t(e')) serves to fix a point in the two microtiles. Of importance is now the 
vector of translation ah between the two points of the microtiles. 

Similarily, an edge in Ei gn will describe a pair of (n + l)-th order microtiles in an n- 
th order microtile. By selfsimilarity again, the corresponding translation vector a e € M. d 
between the two (n + l)-th order microtiles will satisfy 

a e = 0- n a h . (26) 
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Figure 5: A microtile pattern associated with a horizontal arrow e g Eig^ih) (the pattern shown 
has the size of a single tile). The arrow represents the vector r e . 



if e € Eig )Tl (h). See Figure [5] for an illustration. 

Remember from Remark I4.5| that we can identify n oo (^) = II_ 00i o(^)- And the 
inverse of the Robinson map TZ also induces a dense map IP ^ — > t v which is one- 
to-one on the pre- image of 111^ (Cy); we still denote this map by TZ . Hence the 
approximation graph for II" ^ is also an approximation graph for t v . Let Ef denote 
the set of edges whose corresponding paths pass through v at level 0. We may thus adapt 
the above spectral triple to get the spectral triple {C{t v ),9)J , ^lg) (which depends on 
pig) with representation and Dirac defined as in equations (|lip and (fl~2|) by: 

% = P(E? g ) , n 9 (fMe) = f{K-\s(e)))<p(e) , D lg <p(e) = P J g n ^) , e G E? g>n . 

That the bounded commutator axiom is satisfied follows from the following Lemma and 
the fact that Holder continuous functions are dense in C(t v ). 

Lemma 4.6. If f E C(t v ) is Holder continuous (with respect to the Euclidean metric d) 
with exponent a = |°|^ g ^ , then [Di g ,iTi g (f)] is bounded. 

Proof. Suppose / E C(t v ) is Holder continuous with exponent a = r-frar^ > that is, 



d{x,y) 



< C for some C > and all x,y E t v then 



[A ""(/)] II = sup sup 



f{n-\r(e)))-f(1l-\ S {e))) 



d{K-\r{e)),K-\s{e))Y 



d{n- l {r{e)),n- l {s{e))) c 

Pill 



this expression is finite since the first factor is bounded by C. By self-similarity there 
exists C > such that d(TZ~ 1 (r(e)),TZ~ 1 (s(e))) < C'6~ n . And a has been chosen so 
that 0" na p~ n = 1. □ 



We refer to this spectral triple (C(t v ),^ g ,Di g ) as the longitudinal spectral triple ft 



or 



the prototile t v . It should be noted that, although the map TZ 1 is continuous, the 
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topology of t v and rF ^ are quite different and so the Connes distance of this spectral 
triple does not induce the topology of t v . By Theorems 13.41 and 13,81 the longitudinal 
spectral triple has metric dimension 

Slg = J^ff } for all v, but what depends on v is 
the residue of the zeta function. In fact, as compared to the zeta function of the full 
triple it has to be rescaled: Q r = {L v /^2 u L u )Ctr- 

The spectral measure \J( g is easily seen to be the normalized Lebesgue measure on t v , 
as the groupoid of tail equivalence acts by partial translations. 

Similarily to the transversal case we are interested in the quadratic form defined 
formally by 

Ql g (f,9) = V g ([D lg ^ lg {f)]*[D lg ^ lg (g)]) 

v Tr % {\D lg r[D lg ,n lg (f)r[D lg ,n lg (g)^ 
= lml Tv-T-, • ( 27 

sUlg (lg{s) 

Again, this expression has yet little meaning, as we have not yet specified a domain for 
this form. 



4.5 The spectral triple for f2 CT 

We now combine the above triples to get a spectral triple (C(Tl a ),$),D) for the whole 
tiling space The graphs Q and Q have the same set of vertices V, so we notice that 
the identification 

n_oo j+00 (g) = |J n^oCfiO x U-l +oc (g) = (J K,(g) x u^g) 

vev vev 

suggests to construct the triple for Q CT as a direct sum of tensor product spectral triples 
related to the transversal and the longitudinal parts. In fact, IF ^ ((/) x ilg +00 (^) is 
dense in t v x 3t v (see Remark 14. 5 P and so we can use the tensor product construction for 
spectral triples to obtain a spectral triple for C(t v x H^) = C(t v ) <8> C(H^) from the two 
spectral triples considered above. Furthermore, the C*-algebra C(Q, a ) is a subalgebra of 
QXgy C {t v x ^t v ) and so the direct sum of the tensor product spectral triples for the 
different tiles t v provides us with a spectral triple for C(£l): 

£ = 0^r®£f g , vr = 0< r ®^, D = @(Dl®l + X ® D? g ), (28) 
vev vev vev 

where \ ls the grading of the transversal triple (which flips the orientations in Ef r ). The 
representation of a function / € C(£l) then reads 

*(/) = E ftr ® fig . with ftr = <(f) G C(~ tv ) , ff g = G C(t v ) . (29) 

From the results in Section [231 we now get all the spectral information of (C(£l a ), Sj, D). 
To formulate our results more concisely let us call A € non resonant, if A n n ^t°° 

CAe mif for some ca > and non resonant ip £ (0, 2ir) (Definition 13. lip . 
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Theorem 4.7. The above is a spectral triple for C(f2 CT ). Its spectral dimension is 

d log 6 d log 9 

SO = S tr + S ig = h 



log p tr ~ log pig 

and its zeta function Q{z) has a simple pole at sq with strictly positive residue. Moreover, 
for A = ® v A^ r (g> A™ , with either both A\ r and A^ g strongly regular, or one is strongly 
regular and the other the sum of a strongly regular and a non resonant part, then we have 

T{A) = freq^vol^) 7^)7^) . (30) 

In particular, the spectral measure is the unique invariant ergodic probability measure p 
on Q„. 



Proof. Consider first the triple for the matchbox t v x Ht„, that is, the tensor product 
spectral triple for C(£y)<g>C(H^). Applying Lemma[22]we obtain the value Sq = zr^r t + 
- log pi ^ or ^ e aDsc i ssa °f convergence of its zeta function ( v . In particular, this value 
does not depend on v. Furthermore, 

Eoo p / d log 9+2nik \ -p i d log 9—2nik \ 
y . , , r , ■ ^-oo 1 \ -2\og{p tr ~) 1 V -21og(p, g ) ) 

jm(s WW- EuVQl{tu) 2 r(f)log( Ptr )log( Wff ) 

The above number is in fact a strictly positive real number, as it is up to a positive 
factor the mean of two strictly positive periodic functions. It follows that the abscissa of 
convergence for the zeta function of the direct sum of the above triples £ is equal to the 
common value so = Sq. From this we can now determine with the help of Lemma 12.71 
and (|7|) the spectral state. 

HA v tr and A v lg are strongly regular then, by Corollary E2J T{A v tT ®A v lg ) = n v T£(A v tr )T? g (A v lg ) 
with the factor n v = freq(t v )vol(t v ) because the states are normalized. If, say, A\ r is reg- 
ular and A v lg = AJ g reg + Af g nres is the sum of a strongly regular and a non resonant part, 
then 

T(At r (g) Ai g ) = T{A v tr ® Af g sreg ) + T\A\ r <8> Af g nres ) 

lg,nres) 

= n v7t 1 r {A v tT )T2{A v lgsreg ) + 
= n v T^ r (A v tr )Tj" g {A v lg 

,sreg ) + n vT^ r (A v tr )Tj" g (A v lg 

= n v7u.(A^ r )Ti V g (A^ g ) , 

where the second line follows by Corollary 12.81 the third by Lemma 13.101 and the forth 
by Lemma 13.91 (the state of a non resonant operator vanishes: Tj" g {Aig^ nres ) = 0). The 
argument is the same if AJ g is strongly regular, and A\ r the sum of a strongly regular 
and a non resonant part. Hence we get in both cases 

T(A) = t(J2 K ® ai) = n v iz(K)VMD ■ 

V V 

Since n v = freq(i„)vol(i„) is the /i-measure of the matchbox t v x St v we see that the 
spectral measure coincides with \i. □ 
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4.6 Pisot substitutions 



Recall that we consider here aperiodic primitive FLC substitutions which admit a fixed 
point tiling. Their dilation factor 9 is necessarily an algebraic integer. There is a di- 
chotomie: either the dynamical system defined by the tilings is weakly mixing (which 
means that there are no non-trivial eigenvalues of the translation action) or the dilation 
factor 9 of the substitution is a Pisot number, that is, an algebraic integer greater than 1 
all of whose Galois conjugates have modulus strictly smaller than 1. In the second case 
the substitution is called a Pisot-substitution. Let us recall here some of the relevant 
results. We denote by ||x[| the distance of x G R to the closest integer. A proof of the 
following theorem can be found in [3] [Chap. VIII, Thm. I]. 

Theorem 4.8 (Pisot, Vijayaraghavan). Let 9 > 1 be a real algebraic number and a/0 
another real such that \\aO n \\ — > 0. Then we have 

• 9 is a Pisot number. We denote {8j : 1 < j < J} the conjugates of 9 with 9\ = 9 
and \9j + \\ < \9j\. 

• a G Q[0], i.e. a = p a {9) for some polynomial p a with rational coefficients. 

• There exists no such that for n > no 

J 

i=i 

If 9 is unimodular then no = 0. 

We note that p a {9j) ^ for all j, since otherwise p a would be divisible by the minimal 
polynomial of 9 implying that also p a {9) = 0. 

We assume throughout this work that 9 is irrational so that there is at least one other 
conjugate ( J > 1). Note that # _1 is a polynomial (of degree J — 1) in 9 with coefficients 
in jjZ where iV is the constant term in the minimal polynomial for 9. 9 is unimodular 
precisely if N = ±1. 

Recall that a dynamical eigenfunction to eigenvalue (3 £~M. d , is a measurable function 
fp : f2 CT — > C which satisfies fp(oj + t) = e 2n/3 ^ fp(w) for almost all uj G (with respect 
to the ergodic measure p) and all t G M. d . If /« can be chosen continuous then /3 is also 
called a continuous eigenvalue. The set of eigenvalues forms a group which we denote E. 
We call a vector r G M. d a return vector (to tiles) if it is a vector between the punctures 
of two tiles of the same type in some tiling T G £l a . 

Theorem 4.9 (|21j). Consider a substitution with dilation factor 9. The following are 
equivalent. 

1. (3 is an eigenvalue of the translation action. 

2. (3 is an eigenvalue of the translation action with continuous eigenfunction. 
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3. For all return vectors r one has \\(3(r)9 n \\ — > 0. 

In particular, we may assume that all eigenfunctions are continuous, and if there are 
non-zero eigenvalues then 9 must be a Pisot number and for all return vectors (to tiles) 
r and large enough n we have 

J 

exp(2vrz/3(r)r) = exp(-2 m ^^ (r) (%)^) (31) 

i=2 

We would like to give a geometric interpretation of the values of pp^{9j). To better 
illustrate this we consider first only the unimodular situation and explain what changes 
in the general case later. The following theorem can be found in [2|. 

Theorem 4.10. Consider a d- dimensional substitution as above with dilation factor 9. 
If 9 is a unimodular Pisot number of degree J then E is a dense subgroup ofM. d of rank 
dJ. 

Recall for instance from [2j that the maximal equicontinuous factor of (rio-,R rf ) can 
be identified with E the Pontrayagin dual of the group of eigenvalues E equipped with 
the discrete topology. Its M. d action is induced by translation of eigenfunctions, namely 
a r : E ^ E, r eR d , acts as (a r (x))(/3) = e 2m ^xW), X G E, j3 € E. The factor map 
7r : SI — ?• is dual to the embedding of the eigenfunctions in C($l a ). We may choose an 
element T* € CT and consider ir(T*) as the neutral element in E. 

E is a free abelian group of rank d J. We can therefore identify it with a regular lattice 
in an c/J-dimensional Euclidean space W equipped with a scalar product (•,•), and E 
with W/E rec via the map W/E rec 9^h> e 2m ^ 6 E (E rec is the reciprocal lattice). 
With this description of E the dual pairing E x E — > C becomes W/E rec X E 3 (£, f3) (->• 
e 2m \£>P) e C. Furthermore e 27 "'^'^, /3 £ E, is an eigenfunction of the R dJ -action by 
translation on W/E rec to the eigenvalue (-,/3) € M. dJ . Finally, the action of a r becomes 
K(£))C#) = e 2mf3{T ^{P) and hence a r (0) is the unique element of W/E rec satisfying 
/3(r) — (a r (0),/3) € Z for all /3 € E. Now r i— >■ a r (0) is continuous and locally free. It 
follows that the map n->a r (0) lifts to a linear embedding of M d into the universal cover 
W of E. We denote the lift of a r (0) with r. The vector r is thus defined by 

P(r) = V/3e£. 

The image of the embedding, which we denote by U, is simply the lift of the orbit of 
7r(T*). The acting group M d , or equivalently the orbit of 7r(T*), can therefore be identified 
with a subspace U of a euclidean space W. Our next aim is to construct a complimentary 
subspace S to U. 

The endomorphism on M d which is dual to the linear endomorphism 91 on M. d pre- 
serves the group of eigenvalues E cM. d . It therefore restricts to a group endomorphism 
of E. We denote this restriction by (p. With respect to a basis of E, ip is thus an integer 
d J x d J matrix. Denote by </?r its linear extension to W and by ^ the transposqj of 

6 Although taking the transpose is a equalization we have not returned to the original map 81 since 
the two dualizations are with respect to two different dual pairings. 
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<p&. Then we have 

H(r),/3> = <r,w(j3)) = v(P)(r) = OP(r) = 0<r,/3) 

showing that (p^(r) = Or. We know also from [2 J that (p^ has eigenvalues 9j, 1 < j < J 
each with multiplicity J. It follows that U is the eigenspace of <p^ to the eigenvalue 9 
and so by the Pisot-property is the full unstable subspace of ip^. We let S be the stable 
subspace of c^. 

Before we continue the description we provide an example, the Fibonacci substitution 
tiling. Here 9 = is the golden mean and E is the rank 2 subgroup Z[#] C R*. We 

choose the basis {9, 1} for E. Then ip has matrix elements f j ) ^ ^ aS ' °^ course ' 

the same matrix expression. It follows that U is the subspace of W = R 2 generated by 
the vector (9, 1) and S 1 is the subspace generated by (— 9 , 1). It is no coincidence that 
this is the cut & project setup for the Fibonacci tiling, except that we do not have a 
window. 

We return to the geometric description of the polynomicals. Let r be a return vector. 
We choose a lattice base {/3v}v=l— dj f°r E and let {f3 v } v =\... dJ be the dual base for 
E rec . We can express r in the dual base. Since (r, f3 v ) = (3 v {r) we have 

dj 

By Theorem 14.81 there exist polynomials p u with rational coefficients such that 

Pu(r)=p v (9). 

In particular, the vector with coefficients (pv{@)) is an eigenvector of <pL to eigenvalue 
9. The eigenvalue equation Yl^i^Pufi ~ &fivfj,)Pp{0) is a set of dj polynomial equations 
with rational coefficients and hence it is satisfied for 9 whenever it is satisfied for all 
conjugates 9j. Thus (Pu(9j)) is an eigenvector of ip^ to eigenvalue 9j. If j > 1 this vector 
lies thus in the stable manifold S of (pL at 0. We may thus define the following map 
* : {r : r is a return vector} C U — > S, 

J dJ 

j=2 u=l 

The map fi^f* is actually Moody's star map. Indeed, W decomposes into the direct 
sum of vector spaces W = S © U and contains E rec as regular lattice. Let nu and tts 
be the projection onto U and S, respectively, with kernel S and U, respectively. Recall 
that ^2j = iPu(^j) £ Z (we assumed that 9 is unimodular). This can be reinterpreted as 
r + r* G E rec . Since 5 intersects E rec only in the origin^, r* is the unique vector w € S 1 

7 i/5* preserves the intersection £; , " ec n S which is hence invariant under a strictly contracting map, and 
since the intersection is uniformly discrete it must be {0}. 
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which satisfies r + w G E rec . Stated differently, r* = 71$ ° vr^ (r), is the image under 7Ts 
of a preimage in £ ,rec of r under ttu ■ 

Let S2 be the subspace of S generated by the subleading conjugates of 8, i.e. the 
eigenspace of <p^ to the eigenvalues 82, • • • ,$l which are characterised by the property 
that they all have the same modulus: \8j\ = \0z\ for all 2 < j < L < J. We define the 
reduced star map S2 by 

L dj 

3=2 u=l 

By linearity we get 

L 

T,PP(r)(0j) = (^,fi) (34) 
i=2 

which yields the geometric interpretation of the polynomials Ppt r ) we looked for. 

Remark 4.11. If 8 is not unimodular the arguments are principally the same except for 
having to work with inverse limits. The basic differences are that no might be strictly 
larger than in Theorem 14.81 and that <p is no longer invertible. In fact, in the non uni- 
modular case Theorem 14.101 has to be modified in the following way [2]: there exists a 
dense rank dJ subgroup F such that E = lim_ 5 .(i* 1 , (p). In particular, E is an inverse limit 
of dJ-toh. Now one can construct the euclidean space W with its stable and unstable 
subspaces under tp^ as above but for F instead of E. Then E corresponds to the subgroup 
E rec := U n >o n (F rec ) C With these re-interpretations, the map r <— > r* is the 

same, namely r is the lift of a r (0), where a r the action of r £ M rf on the torus F, and r* 
is the unique vector w G S which satisfies r + w G E rec . 

4.7 Dirichlet forms 

The goal of this section is to investigate when the formal expressions for the Dirichlet 
forms can be made rigourous. We will see that, apart from trivial cases, this fixes the 
values for pt r and pi g . The study is technical. We explain the steps of the derivation 
which are similar for both forms, but we do not give all the straightforward but tedious 
details of the technical estimates. 

We start by specializing the formal expression of the Dirichlet form given in equa- 
tions (|2ip and (|22p to our set-up. According to equation (|29p . we view the representation 
of h G C(0 CT ) as the sum of elementary functions ir(h) = X^^eV ^tr ® ^i g . And by The- 
orem [2121 any such h v a is a strongly regular operator on ¥) v a . Hence, by Lemma |2.1U[ we 
can decompose the form as follows 

Q(f,g) = T([D,n(f)]*[D,ir(g)]) = Qi g (f,g) + QMg) , 

8 Note that E rec is defined by this union; it is not the reciprocal lattice of E which wouldn't make 
sense, as E is not a regular lattice in W . 



41 



with 

QM9) = E fre ^)r(l [D? g J? g ]*[Dl,g? g ])lZ(f? r *gl)) (35a) 
vev 

QUf,g) = Y,™i(tv)T(lDiJi]*[Dl,gl] ® i)7S(/5r^) (35b) 

We are going to show in the next two paragraphs, that for /, g in a suitable core, the 
operators [Df , fi g ]*[Df ,gf ] are strongly regular, and the operators [Df r , f£]* [Df r , gf r ] 
are the sums of a strongly regular and a non resonant part. Hence by Corollary 12.81 and 
Lemma 13.101 we get the following decomposition of the forms: 

W, S )^»,«,*)/. (36a) 

Qtr(f,9)=J2 n vQ V M,9tr) [ « <K , (36b) 
vev tv 

where n v = vol(t„)freq(^) is the normalization factor in equation (|30p . d[i v tr is the nor- 
malized invariant measure on (by the results of Section I4.3[) , and d(if g the normalized 
Lebesgue measure on t v (by the results of Section T4.4p . 

Each of the forms Qf and Q" tr are of the type given in Section 13.41 we only have to 
substitute the expression 

§af = /°7e- 1 (r(e))-/o^- 1 (g(e)) _ (3?) 

Pa 

in q" an in equation (|22p for Q^, a G {tr, Ig} and v G V. It follows from Prop. 13.131 that 
the forms Q v tr and Qf g are symmetric, positive definite, and Markovian on the domains 
VI = {/ G Ll(E tv ,dtf r ) : QUfJ) < +00} and Vf g = {/ G : Qf g (f,f) < 

+00}, respectively. 



The longitudinal form Let us first look at the longitudinal part Qf g , which is simpler. 
We show that q" a n (f, g) has a limit for suitable /, g. This will prove that [Df g , /^]* [Df , gf] 
is strongly regular by Corollary 12.51 and will imply the decomposition of Qi g given in 
equation (|36a|) . 

We wish to adapt the parameter p\ g so as to obtain a non-trivial form with core which 
is dense in C l {t v ). 

Given an edge h G Hig, let us denote by ah G M d the translation vector between 
the punctures of the microtiles associated with s(h) and r(h) in the decomposition of 
the tile associated with s 2 (h) = sr(h). If a e G M. d denotes the corresponding vector for 
e G Ef gn (h), equation (|26p gives a e = 9~ n a,h- For n large we can thus approximate 
(uniformly, by uniform continuity of /) 

6 l e 9 f~( K • V)/(x e ) , where x e = Tr x {s{e)) . 

\Pig J 
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Choosing f,g £ C^(t v ), we can substitute the above approximation in qf gn (f,g), up to 
an error term uniform in e, which gives 

//i-i\ 2n i 

where Ej n (h) is the set of edges of type h in E^ g n . 

In order to estimate the above sum, we decompose t v into boxes. Assuming that n 
is large, we choose an integer / = l n such that 1 « I « n, and consider the boxes 
-B 7 = 7£ -1 ([7]) for 7 € ILV(Q), the set of paths of length I which start at v at level 
(here [7] is the set of infinite paths which agree with 7 from level to level The idea 
is that i? 7 should be large enough to allow us to make averages over it, yet small enough 
for continuous functions to be approximately constant on it: F(x e ) ~ i ? (x 7 ), for some 
x 7 € -E> 7 , and all e for which x e £ i? 7 . We have 



Qlg,n(f,9) ~ — EE jfpS ( Qh " V '9)/( X 7) K " Vj ff )s(x 7 ) , 

v p, » J hen lg ^ ^^9," 

for some x 7 € i? 7 , and where Ei gn (h,j) stands for the set of edges of type h whose 
associated (range and source) infinite paths agree with 7 from level to level I. Using 
the estimation of the powers of the substitution matrix in equation ([9]) (note that we use 
the transpose of A for Q here) we have 

#^p,n(^7) L r(7) ^ 2w Ag- 1 -'(l + o(l)) , „, RU1X 

with 

Qg = ( E ^s 2 w) • ( 38 ) 
The last step uses the approximation ji v lg (B 1 )F(x 1 ) ~ f B F(x)d d [i" lg , to obtain 

— q?g,n(f,g) -q v a (f,g) ■■= ci g E fre( i(*s 2 (fc)) / K-v)/ (a h -v)g drf t 

Pig / hc nj, Jtv 



ig 



heH 



Hence, if p/ g = , the sequence qi gn (f,g) converges to q^ g (f,g), therefore by Corol- 
lary [23] we have Q^ g (f,g) = qf g (f,g)- We now have to compute Q\ g from equation (|36a|) . 
summing up the over ti£V, Notice that by the decomposition of the representation 
of a function as vr(/i) = X^eV ^tr ® ^Fg' we have 

E n " / ' ' V )# ^ / ftrgtrdflt r = (a h ■ V ig )f (d h ■ V Jg )p tfyt , 
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where V; ff is the longitudinal gradient on Q a : it takes derivatives along the leaves of the 
folliation; so it reads on the representation of C(Q a ) simply V/ 9 = 1 <8> V K d. Define the 
operator on Ljg(J) CT , dp): 

&lg = cig Vjg/CV^ , with /C = ^ freq(t a a( h ))a h (8)a A . (39) 
We thus we have 

( if, &lg f)i£(n.,dri if Pig = e ~ X 
Ql 9 (f,f) = < H Pig ye' 1 , for all / G Cf^) , (40) 

I +oo if p Jff < 6*- 1 

where C^(r2 CT ) is the space of longitudinally C 2 functions on Q, a . So we see that for 
Pig > is essentially self-adjoint on the domain C 2 g {Sl a ), and therefore the form 

Ql g is closable. For pi g < 6' 1 the form is not closable. 

The transversal form We show that q'^ rn (f,g) decomposes, for suitable /, g, into two 
pieces: one which has a limit, and the other which is oscillating with a phase which will 
be assumed to be non resonant (Definition 13. lip . This will prove that [Dt r , ftr]*[Dtr,9tr] 
is the sum of a strongly regular and a non resonant part, and will imply by Corollary 12,81 
and Lemma 13.91 the decomposition of Qt r claimed in equation (|36bp . 

We wish to adapt the parameter pt r so as to obtain a non-trivial form with core which 
is dense in C(E tv ). 

One might be tempted to consider functions which are transversally locally constant 
to define its core. However, it quickly becomes clear that Qt r vanishes on such functions. 
Indeed, if / 6 C(f2 CT ) is transversally locally constant, then there exists nj G N, such that 
7^ _1 (no jn/ ({?)) gives a partition of H on which / is constant. So we see from equation (|37p 
that <5* r / = for all e G Etr t n, n > n f, and hence Qtr(f, f) = 0. 

We will therefore consider a different core, namely the space generated by the eigen- 
functions of the action. We assume that these form a dense set of L 2 (Cl a ) which is 
equivalent to the fact that the tiling is pure point diffractive. 

Consider an eigenfunction to eigenvalue /3, fp : S1 CT — > C. By definition fa satisfies 
fp{u + r) = e 2n/3 ^ fp(aj) for all cu G £l a and r G R d . Given v G V, h G 7i tr , and any 
e G El n (h), we have ^" 1 (r(e)) = n~ l {s{e))+r e = 7£ _1 (r(e)) + 6 n r h by equation (pi) . 
Then equation (|37p gives 

5 tr fp = {e 2^ m _ l)fpQlz -l {s{e)) 

and hence 

Qtr,n(fl3'fl3') = jl™ 

^tT,n heHtr Ptr 
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As for the longitudinal form, we average over boxes which partition S^: =1Z 1 ([7]), 
for 7 £ Ilo i /(^), and 1 << I << n, and we get 

/ e -27n/3(/>)0™ _ x \ f e 2ml3'(h)8 n _ j\ _ 

Qtr,n(fpJ/3') -Ctr ^ ^ fre q(^2(h)) /_ fpfp'd^r 



heHt 



with 



Ptr 



hen tr 



(41) 



In view of equation (|36bp let us consider 

%r,n(/, 5) = X] "* 1tr,n(ftn9tr) / ffgdlgdtfg ■ 



Summing those terms up over v € V. the integrals yield J n fpfp'dp 
orthogonality of the set of eigenfunctions. Hence we get 



by 



qtr,n(fpi //?')- C tr 5/3/3' II f/3 f ^ freC l(*s 2 (A) ) 



,2m/3{h)e n _ j 



/V 1 



We now use the results of Section f4.6l Using equation (|3ip . expanding the exponential, 

if |0j| < 1 62 1 we approxi- 



and neglecting terms proportional to ^j^rj against 
mate 



e 2 m /3{h)e n _ 1 


2 


n 
Ptr 





i=2 



Ptr 



2n L 



A°>j- a j') n 



ptr j hjl=2 



Here L is such that \0j\ = \02\ for 2 < j < L and we wrote 0j = 1 6*2 1 e*"- 7 ' . The second 
approximation is justified by the fact that Pp(h){9j) 7^ for all j. We assume now that 
the phases are non resonant (Definition 13. 1 1 [) : 

a* - a r + 2<irk + 27r^^fcV > Vk,k'eZ 

fog Pig 

for all j 7^ j'. This implies that the oscillating parts of [Dt r , /tr]* [-Ctn S'ir] are non 
resonant, and hence by Lemma 13.101 do not contribute to T{\JDt r , ftr]*[Dtr, 9tr] ® l) hi 
equation (|35bp . Concerning the non oscillating part, so when j = j', we clearly see that 
they vanish in the limit n — > +00 if pt r > \02\ an d that they tend to +00, as n — > +00, if 
Ptr < \02 1- Finally, if pt r = 1 6*2 1 then the non oscillating part of qtr,n(fp,fp) converges to 



-c tr (27r) 2 ||/ /3 || 2 heq(t s 2 ih) )Y,\Pm(°J 
heHtr 3=2 
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Define the operator At r on the linear space of dynamical eigenfunctions by 

L 

Atrf? = -Ctr(27T) 2 ]T W^(h)) E ^W^i)! 2 & " ( 42 ) 

Then, on the space of dynamical eigenfunctions the transversal form is given by 

{(fp, Atrfp)L 2 (n a ,diJ,) ^ Ptr = \h\ 
i{ Ptr >\9 2 \ . 

+00 if p tr < \9 2 \ 

Clearly, Qt r is closable but trivial if pt r > \O2\1 whereas Qt r is not closable if pt T < \9 2 \. 

Main result We summarize here our results about the Dirichlet forms. For a Pisot 
number 9 of degree J > 1, we denote 9j, j = 2, • • • J, the other Galois conjugates in 
decreasing order of modulus. We write the sub-leading eigenvalues in the form 9j = 
\9 2 \e ia i , 2 < j < L, where ay € [0, 27r). In particular, \9j\ < \9 2 \ for j > L. 

Theorem 4.12. Consider a Pisot- substitution tiling o/R rf with Pisot number 9 of degree 
J > 1. Assume that for all j 7^ j' < L one has 

lo£ \9 I 

ou - a v + 2irk + 2ir ° 1 j f 1 k' / , Vfc.fc'eZ. (43) 
J log 6* 

iSei p/ 5 = # _1 and p tr = |^|- 

If the dynamical spectrum is purely discrete then the set of finite linear combinations 
of dynamical eigenfunctions is a core for Q on which it is closable. Furthermore, Q = 
Qtr + Qi g , and Q tr jig has generator A tr /i g = J2hen tr/lg ^tr/lg 9 iven b V 

tf g fp = -Q 9 (2vr) 2 freq(t s2(h) )/3(a fc ) 2 //3, 
Kfp = -Ctr(27T) 2 freq(t s2(h) ){F h \^ 2 fp 

where c lg = {J2heH lg ^(h))' 1 and c tr = {J2heH tr ^(h))' 1 ■ 

Proof. We have calculated the explicit formulas for and A^ r above. It remains to 
insert (|34|) into (|42p . Eigenfunctions to distinct eigenvalues are orthogonal and hence 
form an orthogonal basis for the Hilbert space by assumption. Standard arguments show 
then that ^tr/lg 1S essen tially self-adjoint and hence the form closable. □ 

Notice that the ratio log(0)/log(|#2|) m equation f|43[) is irrational unless 9 is a Pisot 
number of degree J = 3 and unimodular |22j. 
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4.8 Geometric interpretation of the Laplacians 



We provide an interpretation of the Laplacians as differential operators on the maximal 
equicontinuous factor E of the tiling system. 

Since the dynamical spectrum is pure point and all eigenfunctions continuous the 
factor map ir : £l a — > E becomes an isomorphism between L 2 (f2,/x) and L 2 (E,7]), where 
rj is the normalized Haar measure on E. The Dirichlet form Q can therefore also be 
regarded as a form on L 2 (E,rj). We consider again first the simpler unimodular case in 
which E is a dJ-torus. Each point in E has a tangent space which we may identify with 
U © S, the direct sum of the spaces tangent to the unstable and the stable direction of 
(£>jg, respectively. Now the directional derivative at x along u € U ffi S is given by 

({u,V)fp)(x) = j t fp(x + tu)\ t=0 =2iri{u,(3}fp(x). 

We thus have 

A lg = ci g freq(t s 2 (h) ){a h ,V) 2 , 
At = c tr fre q (t sHh) )(r~ h *,V} 2 . 

In the non-unimodular case one obtains essentially the same: E is an inverse limit of 
tori F = W/F rec with respect to the map ip^, E = lim^(W/F rec ,(p t M ). An element of 
E is a sequence (x n ) n of elements x n € W/F rec satifying x n = <p^(x n +i). The action 
a r of r € M. d on E is thus a r ({x n ) n ) = (x n + ip^ n {r)) n . The continuous functions 
on E axe the direct limit of continuous functions on the tori with respect to the pull 
back of ip l R , C{E) = lim^(C(W/F rec ), The elements of C°°{E) are approximated 

by sequences (f n ) n of elements f n G C°°{W / F rec ) which are eventually modulo the 
equivalence relation identifying (0, • • • , 0, /, — ip^ f, ■ ■ ■ ) with 0. It follows that the 
definition of the directional derivative is compatible with the above equivalence relation 
so that the expression (u, V)(/ n ) n makes sense. In particular, since each eigenfunction 
corresponds to an element of the form (0, • • • ,0, fa, • • • ), f3 £ F, we obtain the same 
formula as above: ((u, V)f/3)(x) = 2m{u, /3)fp(x) and hence also the same formulae for 
Af fl and A|',.. 

This shows that the generators of the Dirichlet forms can be expressed as elliptic 
second order differential operators on E with constant coefficients. 
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